Équation différentielle de Heun et propriétés des fonctions de Heun 


En rédigeant un document concernant la construction de solutions aux problémes extérieures de 
l'équation de Laplace en coordonnées « cyclidiques », je me suis apercu qu'il serait difficile de se 
passer d'une description et de quelques propriétés des fonctions de Heun et de l'équation 
différentielle associée. L'essentiel des résultats donnés ici vient de « NIST Handbook of 
Mathematical Functions : Heun functions ». Je changerais légérement les notations en portant un 
tilde aux traditionnels paramétres a et 6. 


L'équation de Heun s'écrit sous la forme suivante : 
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Contrainte de Fuchs G+B+1=7+8+e 


C'est une équation qui possède 4 points singuliers z=0, z=1, z-a et z-ee. Sur chacun de ces points 
singuliers, si l'on développe une solution en série de Frobenius, l'équation indicielle associée donne 
les exposants et les formes suivantes de développement : 


Développement de Fróbenius 
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Ces développements sont des formes locales valables autour du point singulier dans un rayon de 
convergence défini. Sans entrer dans les détails le rayon de convergence est en général inférieur à 
un. Communément lorsque l'on parle des fonctions de Heun, on parle donc le plus souvent et dans 


un premier temps de ces solutions locales. Elles sont notées : HeunG(a a, à, B y ,ó; z) selon les 


notations de Mathematica, auquel je rajouterais l'indice | pour local, soit HeunG (a q.&, D, y ài) 


Comme ces fonctions sont définies par un développement de Fróbenius, non convergeant partout, 
elles ne sont pas aptes à représenter les solutions d'un probléme aux limites de Sturm-Liouville, par 
exemple entre [0,1] qui sont deux singularités de l'équation de Heun. 


L'équation différentielle de Heun est une équation différentielle Fuchsienne. Tous ces points 
singuliers, y compris le point à l'infini sont des points singuliers réguliers. Les divers paramétres 
sont contraints par la relation de Fuchs qui garanti justement que les points singuliers sont 
réguliers. Pour toute équation fuchsienne du second degré dont les termes peuvent se développer 
autour de la singularité z0 , comme suit : 
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Alors l'équation indicielle est fort simple :" * r(po - 1), 70. 


Pour le point à l'infini, si le développement autour du point singulier à l'infini de P(1/2) et Q(1/C), 
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Appliquons ces résultats à l'équation de Heun : 
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Les développements de Fröbenius sont des développements locaux autour de chaque point 
singulier de l'équation différentielle. De plus la convergence de ces développements n'est assurée 
qu'autour du point singulier dans un rayon qui exclu souvent tout autre point singulier à proximité 
dans le plan complexe. Aussi la résolution d'un problème de type Sturm-Liouville de la construction 
pour une équation différentielle requiert souvent des contraintes supplémentaires portant sur les 
divers paramètres de cette dernière. Lorsque j'évoque le problème de Sturm-Liouville, cela signifie 
classiquement construire la solution de l'équation différentielle dans l'intervalle entre deux points 
singuliers, pour lesquelles sont fixées des conditions limites homogènes portant sur la fonction, sa 
dérivée première ou toute combinaison linéaire de ces dernières. 


Fonctions de Heun analytiques autour d'une, deux ou trois singularités 


En effet aussi bizarre soit-il, l'étude des problémes aux limites de Sturm-Liouville de l'équation de 
Heun n'est pas souvent clairement abordée dans la littérature scientifique, ou alors sous des 
vocables que je n'ai pas bien compris. Par contre les articles débordent de recherches d'autres type 
de développement notamment à l'aide de fonctions hypergéométriques (visant à obtenir soit une 
solution locale autour d'une singularité, soit entre deux singularités, développer dans ce document 
par la suite). C'est un petit paragraphe (31.4) dans le « NIST Handbook of Mathematical 
Functions » qui m'a mis la puce dà l'oreille : « Solutions analytic at Two Singularities : Heun 
Functions ». Il y est trés succinctement décrit la méthode de construction de fonctions analytiques 
sur deux points singuliers, comme 0 et 1, soit peut-être des fonctions aptes à être solutions d'un 
probléme de Sturm-Liouville sur [0,1] aussi simple que celui que je viens de décrire. Toutefois la 
construction proposée s'attache à la résolution d'une équation transcendantale avec une fraction 
continue que je ne suis pas arrivé à implémenter de maniére satisfaisante, car cela divergeait tout 
le temps ! (il s'agit de ma part d'une incompétence temporaire à implémenter convenablement 
l'algorithme de Bouwkanp sur les récurrences à trois termes, mais je me suis rattrapé gráce à 
l'implémentation des fonctions d'ondes sphéroidales réalisé par P.E.Falloon). 


Il existe une dernière sorte de fonctions de Heun, de forme polynomiale, analytique sur les trois 
points singuliers 0,1 et a, il s'agit des polynómes de Heun. 


En réalité la construction de ces trois types de solutions part d'un méme point de départ, la 
constitution de la récurrence du développement de Fróbenius autour de z-0. Pour obtenir les 
développement formels similaires, des secondes solutions on fait appel à des transformations dítes 
F-Homotopiques (un bien grand mot pour un produit de fonction), et pour obtenir les 
développements autour des autres points singuliers de transformations homographiques de 
l'argument z. Mais c'est de tout cela dont nous allons parler maintenant : 


Récurrence du développement de la première solution autour de z=0 
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Le calcul des dérivées première et seconde est alors simple : 
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Pour développer la premiére solution sur les autres points singuliers de l'équation de Heun z-1,a ou 
eo, on utilise des transformations homographiques. 


Pour développer la seconde solution en z=0, on utilise des transformations F-Homotopiques que 
l'on va décrire. 


Pour développer les secondes solutions sur les autres points singuliers de l'équation de Heun z=1,a 
ou ee, on utilise les combinaisons des transformations homographique et F-Homotopiques. 


Transformations F-Homotopiques des solutions de l'équation de Heun 


Il s'agit de rechercher la solution de l'équation de Heun sous trois formes différentes : y(z)-z^'y,(z), 
y(z)-(z-1)?y,(z) et y(z)-(z-a)*yu(z). On démontre alors que les trois fonctions sont également 
solutions de l'équation de Heun mais avec des jeux de paramétres différent : 
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Ces trois transformations homotopiques ont pour propriété simple d'annuler la fonction de Heun 
en chacun des points singuliers z=0,1,a tant que les paramètres y, ô et £ sont des réels inférieurs à 
1: 
»()=2 "ne )= »(0)=0 si y«1 
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Ces trois transformations homotopiques (plus la transformation identité) se combinent deux à 
deux, commutent entre-elles et formes donc un ensemble de 8 transformations de paramètres, 
pas plus de 8 car ce sont des involutions. En appliquant la première transformation on arrive ainsi à 
décrire le développement de la seconde solution autour de z=0 à l'aide d'une fonction de Heun elle- 
même : 
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Les quatre autres transformations homotopiques résultent de combinaisons de celles décrites 
auparavant, dont : 
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Dans le cas d'une formulation à l'aide de variable jacobienne dans l'argument, notamment pour les 
problémes aux limites en coordonnées cyclidiques, les produits qui précédent la fonction de Heun 
ont une valeur plus simple : 
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Pour construire les solutions autour des autres singularités on peut avoir recours aux 
transformation homographique de l'argument z qui permette de permuter des points entre eux. 
Nous allons en décrire simplement deux très simples, mais il en existe 24 en tout, dont on ne va se 
contenter que d'en décrire simplement deux dans ce qui suit. 


Transformations homographiques élémentaires de l'équation de Heun 


La transformation homographique élémentaire z->1-z donne une solution de l'équation de Heun 
suivante, autour de z=1 : 
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Si on lui applique la première transformation F-Homotopique, on obtient la deuxième solution 
autour de z=1, soit : 
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La transformation homographique élémentaire z-»(a-z)/(a-1) donne une solution de l'équation de 
Heun suivante, autour de z-a : 
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Et toujours en appliquant la première transformation F-Homotopique, on obtient la deuxième 
solution autour de z=a, soit : 
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Ce sont les 6 formes de solutions données dans le « NIST Handbook of Mathematical Functions » , 
formules 31.3.5 à 31.3.9. 


Dans une article de R.S.Maier de 2006, « The 192 Solutions of the Heun Equation », l'auteur donne 
le tableau exact de toutes les combinaisons possibles entre les 24 types de transformations 
homographiques et les 8 types de transformations F-Homotopiques pour construire les solutions 
locales de l'équation de Heun, aux 4 points singuliers. Cela donne 24x8=192 solutions possibles. La 
table 2 de l'article les énumère. Elles sont d'une grande importance y compris pour la construction 
des solutions bi-analytiques en deux points singuliers (valable pour les problèmes de Sturm- 
Liouville) car on peut toujours appliquer ces transformations pour atteindre d'autres solutions. 


Récurrence du développement de la premiére solution autour de z-1 


Il est également intéressant de réaliser le même développement autour de la seconde singularité 
z-1. Cela permet d'illustrer l'effet des transformations homographiques. 
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Le calcul des dérivées premiére et seconde est alors simple : 
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On retrouve bien l'effet de la transformation homographique z-»1-z En appliquant cette 
transformation des paramétres comme suit : 
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sur la récurrence initiale du développement autour de z=0, on retrouve la récurrence du 
développement autour de z-1. 


Récurrence du développement de la premiére solution autour de z-a 


Réalisons le méme développement autour de la seconde singularité z-a. Cela permet d'illustrer 
l'effet des transformations homographiques. 
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Première solution y, (z) = HeunG, (a. q.à, f, 7,0; z)= D c; 
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Le calcul des dérivées premiére et seconde est alors simple : 
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On retrouve bien l'effet de la transformation homographique z-»(a-z)/(a-1). En appliquant cette 
transformation des paramétres comme suit : 
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sur la récurrence initiale du développement autour de z-0, on retrouve la récurrence du 
développement autour de z-a. 


Solutions analytiques aux deux points singuliers z=0 et z-1, récurrence du développement, 
valeurs propres matricielles et condition d'existence de la solution analytique en z-0 et z-1 


La construction que je propose est évoquée vaguement dans la littérature, mais à chaque fois 
éludée lorsqu'il s'agit d'une description concréte. Aussi je vais tenter d'exprimer la récurrence et la 
détermination des coefficients du développement comme un probléme aux valeurs propres d'une 
matrice infinie. Revenons au développement autour de z=0, il peut s'exprimer comme un système 
d'équation linéaire de la forme : 
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Autrement dit, il s'agit de trouver les valeurs propres q de la matrice M, et une fonction analytique 
en deux points posséde un développement non trivial et convergeant si q est une valeur propre de 
cette matrice M (soit également que le déterminant de la matrice M-ql s'annule). Les vecteurs 
propres de la matrice, rangés par ordre croissant des valeurs propres contiennent chacun les 
coefficients du développement approchée de la fonction bi-analytique aux deux points z=0 et z=1. 
Plus la dimension de la matrice est importante, meilleure est la convergence des premiéres valeurs 
propres. Pour ma part j'ai préféré ce type de construction à celles proposées dans le paragraphe 
31.4 consistant à résoudre une équation transcendantale impliquant la fraction continue que l'on 
peut également construire avec la récurrence précédente. 


Schéma de construction approchée des fonctions de Heun Bi-analytiques en deux points 
singuliers de l'équation 


On a donc choisi systématiquement la forme algébrique de l'équation de Heun. Pour une fonction 
bi-analytique de Heun, voici le schéma de construction des valeurs propres et fonctions propres 
pour le probléme de Sturm-Liouville entre [0,1] telles que la fonction et sa dérivée premiére soient 
finies aux limites O et 1 (c'est à dire analytique en ces deux points et avec Mathematica, la 
résolution numérique est construite comme suit : 
n1-20; 
HeunEquationNuAlgebriques[a ,a ,6 ,y ,8 ,x ]:= 

X (x-1) (x — a) Derivative[2][y][x]* x (x-1) (x-a) (y/x *6/(x-1)*( a +6+1-y -6)/(x-a)) Derivative[1][y][x]*x a 6 y[x]; 
{HeunEigenValuesAlgebriques, HeunEigenFunctionsAlgebriques} = 

NDEigensystem[ HeunEquationNuAlgebriques[a, a, 6, y, 6,x], y[x], (x, 0, 1}, n1, 

Method -> ("PDEDiscretization" -> ("FiniteElement", ("MeshOptions"-» ("MaxCellMeasure" -> 0.0001}}}}]; 
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Pour la résolution numérique avec Mathematica et sa comparaison avec la construction 
matricielle, voici un exemple de graphe obtenue pour la fonction de Heun avec les paramétres 
suivant : a = 1/0.8; a= 1/3; 6 = 1; y= 2/3; 8 = 1; £ =a+6+1-y-6; la 4ème valeur propre : 


HeunEigenFunctionsAlgebriquesIn] 
f1 


LesFonctionsPropresIn] 
g1 


-0.2 


-0.4 


Et les deux jeux de valeurs propres obtenues numériquement et par la méthode matricielle : 


(0.157896,-0.860301,-3.07963,-6.50981,-11.1508,-17.0025,-24.0649,-32.338,-41.8219,-52.5164,-64.4217,-77.5376,- 
91.8643,-107.402,-124.15,-142.108,-161.278,-181.658,-203.249,-226.051} 
(0.157896, -0.860303,-3.07967,-6.51018,-11.1536,-17.0191,-24.1413,-32.6143,-42.618,-54.3972,-68.2073,-84.2813,- 
102.829,-124.051,-148.145,-175.314,-205.766,-239.72,-277.404,-319.053} 


Remarque : il est clair que les dérivées premières de la fonction ne s'annulent pas aux bords de 
l'intervalle, donc ce n'est pas à proprement parler un problème de Neumann. 


Pour résoudre un problème de Sturm-Liouville avec des conditions aux limites homogènes de 
Dirichlet ou Dirichlet-Neumann sur l'un des bords du domaine, il suffit d'appliquer les 
transformations F-Homotopique Lll ou Ill, ou leur combinaison pour prendre en compte une 
solution qui s'annule nécessairement à l'une des bornes de l'intervalle [0,1], avec la contrainte des 
paramètres y, ô et & étant des réels inférieurs à 1. Cela veut concrètement dire que l'on va 
construire la matrice M avec des paramétres de Heun modifiés par la transformation F- 
Homotopique correspondante. 


Si l'on travaille sur un intervalle quelconque [51,52] où s1, s2 est l'une quelconque des singularités 
parmi {0,1,a,}, alors il suffit (mais avec soin) de combiner également avec une transformation 
homographique de l'argument z. Toutes les combinaisons possibles de ces transformations sont 
remarquablement résumées dans la table 2 de l'article « The 192 Solutions of the Heun Equation » 
dans laquelle figure les expressions explicites des nouveaux paramétres de Heun. 


Avec les Notation NIST modifiée hi(y,8, 5; z)HeunG,, aq, & Br 0: (z)) n —1,2,...,007 


Lorsque que l'on applique une transformation homographique à l'argument, tous les paramètres 


de Heun sont susceptibles d'étre modifiés. 


Pratiquement voyons ce que cela donne pour les divers types de problémes aux limites que l'on 
peut rencontrer : 


Intervalle z€[0,1]: Prenons le cas où la fonction bi-analytique doit s'annuler en z=0 et être 
analytique en z-1, alors on applique l'algorithme avec les substitutions suivantes de paramétres : 


Algorithme (q, y(2)) = HeunBiAnalytiquela, +1-—7, B -1—-y,2-y,6, z) 
e=&+ß+1-y-ô 
Valeur propre d -4-(aó +e) y)=q la5+ã -B+1-y si y) 
y(z)- z ylz) 


Pour la résolution numérique avec Mathematica et sa comparaison avec la construction 
matricielle, voici un exemple de graphe obtenue pour la fonction de Heun avec les paramètres 
suivant : a = 1/0.6; a= 1/5; 6 = 1; y= 1/6; 6 = 2/3; £ =æ+6+1-y-6; avec la 4ème valeur propre : 


1.0 


HeunEigenFunctionsAlgebriquesIn] 
f 


LesFonctionsPropresIn] 
g1 


-0.5 


I: 


Le jeux de valeurs propres est le suivant: {-0.63782,-3.31564,-8.18189,-15.2139,-24.4103,- 
35.7707,-49.2951,-64.9835,-82.8359,-102.852,-125.032,-149.376,-175.885,-204.557,-235.394,- 
268.401,-303.592,-341.013,-380.786,-423.156,-468.513,-517.332,-570.071,-627.102,-688.719,- 
755.161,-826.648,-903.394,-985.62,-1073.56,-1167.45,-1267.54,-1374.1,-1487.42,-1607.78,- 
1735.49,-1870.88,-2014.27,-2166.03,-2326.53,-2496.15,-2675.32,-2864.46,-3064.03,-3274.53,- 
3496.45,-3730.34,-3976.78,-4236.37,-4509.78,-4797.69,-5100.86,-5420.07,-5756.18,-6110.12,- 
6482.88,-6875.55,-7289.29,-7725.39,-8185.25,-8670.41,-9182.57,-9723.62,-10295.7,-10901.1,- 
11542.6,-12223.1,-12946.3,-13716.1,-14537.3,-15415.6,-16357.9,-17372.5,-18470.1,-19664.2,- 
20973.3,-22422.9,-24051.2,-25920.7,-28150.1] et la construction par le développement de type 
Fróbenius respecte sur tout l'intervalle l'équation différentielle de Heun. 


Intervalle z€[0,1] : Lorsque la fonction bi-analytique doit s'annuler en z-1 et étre analytique en 
z=0, alors on applique l'algorithme avec les substitutions de paramètres suivantes : 


Algorithme (4,»(z))= HeunBiAnalytiquela, à +] 8,B +1-0,7,2 5,2) 
s=à+B+1-y-6 
Valeur propre d 2 q-ay(1-6) 


3()--2)"^ yl) 


jeee deem jer 


»(z)- Gor AGE Sioa y," je Y 6-ye;z 
j-0 


j=0 j=0 


X()-(-z y)  ».G)--2)7 v'e)-U-8X1-2)° y.) 


y", ()- 27," )- 20-80 -z)" y, G)« 8-617 z)" », (0) 


j-2 


Pour la résolution numérique avec Mathematica et sa comparaison avec la construction 
matricielle, voici un exemple de graphe obtenue pour la fonction de Heun avec les paramètres 
suivant : a = 1/0.3; a= 1/5; 6 =7/8; y= 2/3; 6 = 5/6; £ =a+6+1-y-6; avec la 3ème valeur propre : 


I 


1.0 


HeunEigenFunctionsAlgebriquesIn] 
0 ft 


LesFonctionsPropresIn] 
g1 


-1.54 


Le jeux de valeurs propres est le suivant: {-0.224463,-5.36739,-16.1004,-32.4335,-54.3665,- 
81.8995,-115.032,-153.765,-198.098,-248.031,-303.564,-364.696,-431.429,-503.761,-581.694,- 
665.226,-754.358,-849.09,-949.423,-1055.35,-1166.89,-1284.02,-1406.75,-1535.08,-1669.01,- 
1808.55,-1953.68,-2104.41,-2260.74,-2422.67,-2590.2,-2763.34,-2942.07,-3126.43,-3316.46,- 
3512.36,-3714.61,-3924.27,-4143.06,-4373.11,-4616.4,-4874.45,-5148.27,-5438.6,-5746.07,- 
6071.29,-6414.92,-6777.67,-7160.36,-7563.83,-7989.05,-8437.04,-8908.92,-9405.9,-9929.31,- 
10480.6,-11061.3,-11673.1,-12317.8,-12997.6,-13714.6,-14471.4,-15270.7,-16115.6,-17009.6,- 
17956.5,-18960.9,-20027.8,-21163.3,-22374.4,-23669.3,-25058.1,-26553.2,-28170.,-29928.8,- 
31856.2,-33990.,-36386.3,-39136.9,-42416.1} et la construction par le développement de type 
Fróbenius respecte sur tout l'intervalle l'équation différentielle de Heun. 


Intervalle z€[0,1] : Lorsque la fonction bi-analytique doit s'annuler en z=0 et z-1, alors on applique 
l'algorithme avec les substitutions de paramétres suivantes : 


Algorithme (q, y(z)) = HeunBiAnalytique(a, & -y-ó- 2,B 7-8+2,2-7,2 5,z) 
e=à+B+1-y-6 
Valeur propre q= e(y = 1)+ aly +ô- 2)+ q 


Pour la résolution numérique avec Mathematica et sa comparaison avec la construction 
matricielle, voici un exemple de graphe obtenue pour la fonction de Heun avec les paramètres 
suivant : a = 1/0.6; a= 1; 6 = 1/5; y= 2/3; 6 = 5/6; £ =a+8+1-y-6; avec la 3ème valeur propre : 


HeunEigenFunctionsAlgebriquesIn] 
f 


LesFonctionsPropresIn] 
g1 


Le jeux de valeurs propres est le suivant: {-0.45801,-3.15801,-8.02651,-15.0592,-24.2558,- 
35.6164,-49.1409,-64.8293,-82.6817,-102.698,-124.878,-149.222,-175.73,-204.402,-235.238,- 
268.238,-303.402,-340.73,-380.221,-421.878,-465.7,-511.697,-559.89,-610.343,-663.204, - 
718.749,-777.386,-839.592,-905.813,-976.401,-1051.62,-1131.66,-1216.71,-1306.92,-1402.46,- 
1503.5,-1610.21,-1722.79,-1841.43,-1966.32,-2097.7,-2235.77,-2380.77,-2532.95,-2692.54,- 
2859.81,-3035.04,-3218.49,-3410.47,-3611.29,-3821.24,-4040.68,-4269.94,-4509.38,-4759.39,- 
5020.34,-5292.67,-5576.78,-5873.15,-6182.24,-6504.55,-6840.61,-7190.97,-7556.22,-7936.97,- 
8333.88,-8747.63,-9178.96,-9628.66,-10097.5,-10586.5,-11096.5,-11628.6,-12183.8,-12763.3,- 
13368.5,-14000.6,-14661.3,-15352.1,-16074.8,-16831.4,-17624.2,-18455.5,-19328.,-20244.9,- 
21209.5,-22225.8,-23298.4,-24432.4,-25634.1,-26910.9,-28271.5,-29726.9,-31291.,-32981.5,- 
34822.3,-36847.,-39105.5,-41680.,-44726.6] et la construction par le développement de type 
Fróbenius respecte sur tout l'intervalle l'équation différentielle de Heun. 


Intervalle z€[a,1] avec 0«a«1 : Prenons le cas où la fonction bi-analytique doit être analytique en 
z=a et z-1, alors on applique l'algorithme avec les substitutions suivantes de paramètres : 


Algorithme (q, y) = HeunBiAnalytique| = a,ã,ßB,ô, yA- z) 
£-üf«1l-y-ó 
Valeur propre ÿ -à B -q 


j=+% j=+% 


net) OE el) OÈ -Denli 


j=0 j=0 


Pour la résolution numérique avec Mathematica et sa comparaison avec la construction 
matricielle, voici un exemple de graphe obtenue pour la fonction de Heun avec les paramètres 
suivant : a = 0.33; a= 1/3; 6 = 1; y= 2/3; 6 = 5/6; € =a+6+1-y-6; avec la 7ème valeur propre : 


HeunEigenFunctionsAlgebriquesIn] 
fi 


LesFonctionsPropresfn] 
g1 


Le jeux de valeurs propres est le suivant : (0.209512, 1.20428, 3.39242, 6.77412, 11.3494, 17.1183, 
24.0819,32.245,41.6286,52.3002,64.4145,78.2232,94.0219,112.088,132.663,155.967,182.212,211 
.62,244.422,280.863,321.204,365.725,414.725,468.525,527.47,591.934,662.318,739.061,822.64,9 
13.574,1012.44,1119.86,1236.54,1363.25,1500.87,1650.38,1812.9,1989.71,2182.3,2392.4,2622.0 
7,2873.79,3150.61,3456.32,3795.84,4175.71,4605.02,5097.33,5674.65,6378.55,7314.37} et la 
construction par le développement de type Fróbenius respecte sur tout l'intervalle l'équation 
différentielle de Heun. 


Intervalle z€[a,1] avec 0«a«1 : Prenons le cas où la fonction bi-analytique doit être analytique en 
z-1 et s'annuler en z-a, alors on applique l'algorithme avec les substitutions suivantes de 
paramétres : 


Algorithme (q, y (z)) = HeunBiAnalytique( a, &+y+ 8,-B +y +0,0,7,1- z) 
e=&+B+1-y-6 
Valeur propre. d -àf +6(1 —e)-q 
y(z)- za) y) 


J=+0 j=+0 J=+0 


xz) Zc,f-2)  p'G)=-Xic,f-2)" e= Y 60-0)e,,-2)7 


j-0 j-0 j= 


5)=6-0)" ze) sG-6G-a» GC» e) le-a) 


dz 


2 


r O=- 7,023, Ele-alea) 


dz 


Pour la résolution numérique avec Mathematica et sa comparaison avec la construction 
matricielle, voici un exemple de graphe obtenue pour la fonction de Heun avec les paramètres 
suivant : a = 0.2; a= 1/3; 6 = 6/7; y= 2/3; 6 = 5/6; £ =æ+6+1-y-6; avec la 7ème valeur propre : 


HeunEigenFunctionsAlgebriquesIn] 
fi 


LesFonctionsPropresfn] 
g1 


-0.5 


Le jeux de valeurs propres est le suivant : {0.262529, 1.3044, 3.3107, 6.28555, 10.229, 15.1412, 
21.0231,27.8788,35.7235,44.6012,54.6066,65.8915,78.6382,93.0224,109.195,127.287,147.417,16 
9.702,194.257,221.202,250.658,282.751,317.611,355.372,396.174,440.159,487.478,538.284,592. 

737,651.006,713.261,779.683,850.46,925.786,1005.87,1090.91,1181.15,1276.8,1378.12,1485.36, 

1598.79,1718.7,1845.37,1979.12,2120.28,2269.21,2426.28,2591.87,2766.4,2950.32,3144.11,3348 
.28,3563.35,3789.94,4028.65,4280.17,4545.25,4824.67,5119.32,5430.14,5758.2,6104.65,6470.79, 
6858.04,7268.03,7702.58,8163.77,8653.98,9176.,9733.08,10329.1,10968.7,11657.6,12403.,13214 
.3,14104.,15089.4,16196.7,17468.4,18985.5,20954.4} et la construction par le développement de 
type Fróbenius respecte sur tout l'intervalle l'équation différentielle de Heun. 


Intervalle z€[a,1] avec 0«a«1 : Prenons le cas où la fonction bi-analytique doit être analytique en 
z=a et s'annuler en z-1, alors on applique l'algorithme avec les substitutions suivantes de 
paramétres : 


Algorithme (q, y (z)) = HeunBiAnalytique( a,à — Ô + LB ô+1,2-6,7,1 z) 


e=&+B+1-y-6 
Valeur propre 3=(&-5+1)8-6+1)+ ay(6 —1)-q 


z 
b 
E] 
—. 
N 
— 
Il 
l 


Sjel- O= E JD, (1-2) * 


j=0 j=0 


zel- 3,0) sG-0-27» C» 6) 1-27) 


dz 


y,(z)- c, dl-z 


` 
Il 
© 


Pour la résolution numérique avec Mathematica et sa comparaison avec la construction 
matricielle, voici un exemple de graphe obtenue pour la fonction de Heun avec les paramètres 
suivant : a = 0.33; a= 1/3; 6 = 1; y= 2/3; 6 = 5/6; € =a+6+1-y-6; avec la 7ème valeur propre : 


HeunEigenFunctionsAlgebriquesIn] 
fi 


LesFonctionsPropresfn] 
g1 


Le jeux de valeurs propres est le suivant : (0.297074, 1.48616, 3.87315, 7.45377, 12.228, 18.1959, 
25.3587,33.7225,43.3124,54.2048,66.5681,80.6635,96.7866,115.21,136.17,159.884,186.565,216. 
434,249.724,286.679,327.562,372.655,422.256,476.688,536.297,601.459,672.577,750.09,834.476 
,926.26,1026.01,1134.37,1252.03,1379.78,1518.48,1669.14,1832.86,2010.94,2204.87,2416.39,26 
47.57,2900.89,3179.41,3486.96,3828.45,4210.44,4642.1,5137.,5717.26,6424.6,7364.82] et la 
construction par le développement de type Fróbenius respecte sur tout l'intervalle l'équation 
différentielle de Heun. 


Intervalle z€[a,1] avec 0«a«1 : Prenons le cas où la fonction bi-analytique doit s'annuler en z=1 et 
z-a, alors on applique l'algorithme avec les substitutions suivantes de paramètres : 


Algorithme (q, y(z)) = HeunBiAnalytiquell —a,-à+y+ 1-8 +7+12-06,7,1- z) 
e=à+p+1-y-6 
Valeur propre q= (à y Jd y 1)+ yle 1)+ay(8 -1)-q 


y(z)- (z-a) = (1-2) xG) 
j=+% Jj j=+% 


ne) 3 e,-zy ne-E la a 3 j0-0e,, 0-2)" 


j-0 j-0 j-0 


Pour la résolution numérique avec Mathematica et sa comparaison avec la construction 
matricielle, voici un exemple de graphe obtenue pour la fonction de Heun avec les paramétres 
suivant : a = 0.33; a= 1/3; 6 = 1; y= 2/3; 6 = 5/6; € =a+6+1-y-6; avec la 5ème valeur propre : 


HeunEigenFunctionsAlgebriquesIn] 
0 " 


LesFonctionsPropresInl 
g1 


Le jeux de valeurs propres est le suivant : (0.408648, 1.80104, 4.38704, 8.16658, 13.1397, 19.3065, 
26.6681,35.23,45.0161,56.1008,68.6514,82.932,99.2446,117.867,139.039,162.979,189.901,220.0 
26,253.585,290.825,332.008,377.416,427.349,482.129,542.103,607.647,679.166,757.099,841.926 
,934.171,1034.41,1143.27,1261.46,1389.77,1529.06,1680.33,1844.7,2023.46,2218.11,2430.38,26 
62.36,2916.54,3195.97,3504.49,3847.03,4230.16,4663.05,5159.32,5741.13,6450.3,7392.84} et la 
construction par le développement de type Fróbenius respecte sur tout l'intervalle l'équation 
différentielle de Heun. 


Intervalle z€[1,a] avec a»1: alors les 4 implémentations de l'algorithme s'obtiennent en 
substituant à chaque les expressions (z-a) et (1-z) en (a-z) et (z-1). 


Les 4 dernières implémentations que je vais traiter sont celles sur un intervalle [0,a] : 


Intervalle z€[0,a] avec a«1 : Prenons le cas où la fonction bi-analytique doit être analytique en z=0 
et z=a, alors on applique l'algorithme sans substitution de paramètres : 


Algorithme (a. y(z)) = HeunBiAnalytique(a, &, B, 76,2) 


s=à+B+1-y-6 
Valeur propre q 


j=+% jt Je 


*u)9»e 5g "ie X-i), 
j-0 


j-0 j-0 


Pour la résolution numérique avec Mathematica et sa comparaison avec la construction 
matricielle, voici un exemple de graphe obtenue pour la fonction de Heun avec les paramétres 
suivant : a = 0.8; a= 1/3; 6 = 1; y= 2/3; 6- 1; € =a+6+1-y-6; avec la 4ème valeur propre : 


1.0 


HeunEigenFunctionsAlgebriquesIn] 
f 


LesFonctionsPropresIn] 
g1 


-0.5 


I: 


Le jeux de valeurs propres est le suivant: {0.144416,-0.503998, -2.11 72,-4.7005,-8.25546,- 
12.7931,-18.3502,-25.0148,-32.9382,-42.3108,-53.3247,-66.1561,-80.9678,-97.9176,-117.164,- 
138.868,-163.198,-190.329,-220.442,-253.727] et la construction par le développement de type 
Fróbenius respecte sur tout l'intervalle l'équation différentielle de Heun. 


Intervalle z€[0,a] avec a«1 : Prenons le cas où la fonction bi-analytique doit s'annuler en z=0 et 
étre analytique en z-a, alors on applique les mémes substitutions de paramétres que pour 
l'intervalle [0,1] avec les mêmes conditions limites soit : 


Algorithme (q, y(2)) = HeunBiAnalytiquela, +1-—7, B -1-y,2-y,6, z) 
e=G+B+1-y-6 
Valeur propre d -q-(a8-cXI- y) q- a6 -à - B--1-y -6f0- y) 
y(z)- z*y(G) 


»()-YXeuz  vx(G)Yez" y le) XG- 
j=0 j-0 j=0 


Pour la résolution numérique avec Mathematica et sa comparaison avec la construction 
matricielle, voici un exemple de graphe obtenue pour la fonction de Heun avec les paramètres 
suivant : a = 0.8; a= 1/3; 6 = 1; y= 2/3; 6 = 1; £ =a+6+1-y-6; avec la 8ème valeur propre : 


4r 


HeunEigenFunctionsAlgebriquesInl 
f 


LesFonctionsPropresfn] 
g1 


Le jeux de valeurs propres est le suivant : (0.0219058,-0.934107,-2.87055,-5.77719,-9.65732,- 
14.5272,-20.4357,-27.4895,-35.8547,-45.7249,-57.2876,-70.7139,-86.1634,-103.793,-123.76,- 
146.227,-171.363,-199.342,-230.348,-264.572,-302.213,-343.481,-388.598,-437.795,-491.316,- 
549.421,-612.382,-680.488,-754.047,-833.386,-918.854,-1010.82,-1109.69,-1215.88,-1329.86,- 
1452.13,-1583.21,-1723.7,-1874.23,-2035.5,-2208.27,-2393.39,-2591.78,-2804.5,-3032.7,- 
3277.72,-3541.04,-3824.4,-4129.8,-4459.56,-4816.46,-5203.85,-5625.84,-6087.57,-6595.69,- 
7159.07,-7790.09,-8507.17,-9340.37,-10346.5) et la construction par le développement de type 
Fróbenius respecte sur tout l'intervalle l'équation différentielle de Heun. 


Intervalle z€[0,a] avec a«1: Lorsque la fonction bi-analytique doit s'annuler en z-a et être 
analytique en z=0, alors on applique l'algorithme avec les substitutions de paramètres suivantes : 


Algorithme (q,y(z))= HeunBianalytiquela, ü-ytó, B Fy4 5,7.6,2) 
s=G+B+1-y-6 
Valeur propre q=q -y(1-«) 


3)-(a- z)* yl) 


j=+0 j=+0 j=+0 


yz)» ON d y,'(z)= ML y," (z)- XU- e,z 


j=0 j=0 j=0 


Pour la résolution numérique avec Mathematica et sa comparaison avec la construction 
matricielle, voici un exemple de graphe obtenue pour la fonction de Heun avec les paramètres 
suivant : a = 0.8; a= 1/3; 6 = 1; y= 2/3; 6 = 1; £ =a+6+1-y-6; avec la 4ème valeur propre : 


L 


HeunEigenFunctionsAlgebriquesIn] 
f 


LesFonctionsPropres{n] 
gi 


Le jeux de valeurs propres est le suivant : (0.050468,-0.933045,-2.87042,-5.77642,-9.65322,- 
14.5099,-20.3793,-27.3459,-35.5629,-45.2299,-56.5493,-69.704,-84.86, -102.176,-121.81,- 
143.926,-168.689,-196.276,-226.867,-260.655} et la construction par le développement de type 
Fróbenius respecte sur tout l'intervalle l'équation différentielle de Heun. 


Intervalle z€[0,a] avec a«1 : Lorsque la fonction bi-analytique doit s'annuler en z=0 et z-a, alors on 
applique l'algorithme avec les substitutions de paramétres suivantes : 


Algorithme (q, y(2)) = HeunBiAnalytiquela.- à +0 + 1-f +0 +1,2-7,6, z) 
e=à+B+1-y-6 
Valeur propre d -&*y -2*aó(y -1)*q 


y()- z*(-:)7»G) 


Pour la résolution numérique avec Mathematica et sa comparaison avec la construction 
matricielle, voici un exemple de graphe obtenue pour la fonction de Heun avec les paramétres 
suivant : a = 0.8; a= 1/3; 6 = 1; y= 2/3; 6 = 1; £ =a+6+1-y-6; avec la 4ème valeur propre : 


HeunEigenFunctionsAlgebriquesIn] 
f 


LesFonctionsPropresIn] 
g1 


Le jeux de valeurs propres est le suivant: {-0.178936,-1.47145,-3.73139,-6.96056,-11.162,- 
16.3494,-22.5675,-29.9202,-38.578,-48.7446,-60.6172,-74.3722,-90.1717,-108.173,-128.534,- 
151.418,-176.992,-205.433,-236.923,-271.654] et la construction par le développement de type 
Fróbenius respecte sur tout l'intervalle l'équation différentielle de Heun. 


Solutions polynomiales analytiques aux trois points singuliers z=0 z-1 et z=a, récurrence finie du 


développement, valeurs propres matricielles 


L'étape suivante des solutions de l'équation est celle pour laquelle la matrice M est de dimension 
finie, car les termes de la récurrence s'annulent au delà d'un certain nombre n. Dés lors on 
construit nécessairement une solution polynomiale de degré n pour laquelle il y a exactement n 
valeurs propres qnm de la matrice M, Pour que la récurrence soit finie une condition 
supplémentaire sur un des deux paramétres a ou 6 qui doit étre un entier négatif. Cela donne un 
schéma de construction de la matrice à diagonaliser suivant : 


R; = a(j+1Y7 + j) 


Q,=j{-1+y7Xa+1)+e+a6} | x 
5. j) Q--n ou p =-n Rc;\-(O,+ak,+Pe,,=0 


Pa — 0> Cua ES 0 


c,=0 oc-1 


-Q R, 0 
P -Q R 0 " 0 
0 Zr m E 0 is 0 
Ce sese M, =| 2: 0 Pas Qm R3 0 ie (M, -t,)C, - 0 
p^ au^ HS Ba cms. Roa 0 
0 0 . 0 Pa Os Rm 
| 0 "E "T 0 P, 2501] 


q,, valeurs propres de M, (a.a. B. y. 5) P5 (z)- Sc, (a? e m € {0,1,2,..,n} 


n,m j,n,m 
j-0 


Notation | valeurs propres VM, (a;à. 8.7.6) et P’ (z) HeunG,(a,q°,.&.f,7,6:z) 


Compte tenu des huit transformations F-Homotopiques laissant invariante l'équation de Heun, on 
peut construire huit solutions apparentées « polynomiales », comme suit. 


Commencons par la premiére transformation F-Homotopique H1 avec le schéma de construction 
suivant : 


q 7 4, -(aó « efl - y) 
H, 4à,-à-4l-y B,=B+1-7 contrainte @,=-n=@+l-y=a@=y-l-n 


V1 =2-Y ô, =ô €; 5E 
P! (z)= z" B^ (z)— q!, eVp|M, (aia, =-n,8,,1,.6,) et P! (z)= HeunG,(a.q! , Gi, B, y,.6,:) 


n,m n,m 


La deuxième transformation F-Homotopique H2 donne le schéma de construction suivant : 


q=q-ay(1-ô) 
H, s&=@+1-ő6 B,-B-1-8 contrainte &,=-n=à&+1-5—=&-6-1-n 


Vn =r Ôr =2- 8 Ep —€ 


P (z)- (T n. (z)= Gin € VM, a.i, = -n, Bn: Y nn ) et [s (z)- HeunG, an 8, Burns 8n32) 


La troisiéme transformation F-Homotopique H3 donne le schéma de construction suivant : 


qum ^ q*Y(0.-6) 
H, 4ügy =&+1-E B,-B-1-s& contrainte Gy=-n=&+1-e—à@=Ee-1-n 
Yum Y Ôn = Ô £g = 2€ 


Et &+ßB+l=y+ő8 +e >8=&+ß+1 y e=f nY 


l-e 
p (z) = í = z) p (z) = Vd € Vp|M, (a: y A, Bu > Var Our ) et ps (z) - HeunG, (a. UE Q > Ba » Vin» Ôm? z) 


Les 4 autres transformations donnent les constructions suivantes de solutions assimilées 
polynomiales : 


Gin =q+ell-y)-aļy+8-2) 

H, 4àj,7ü*2-y-Ó — f,y-*2-y-Ó contrainte à,,2-n-ü*2-y-ó2ü-y*ó-2-n 

Yın=2-Y O2 0 Eye 

Per B) att erp fs ora) et Ee) temas sara 
Qim=q+i-e+(1-yf1+aô) 


Hj, id, -ü42-y-£6 Brun 7 B*2-y-6 contrainte d,,, 2-nzüc2-y-£2üzsyte-2-n 


lui d Ó m =ô &g 72-6 


l-e 
LH) yl Z| DLI LII "IE. BLU \ ED . 
Ts (z)=z ud Po ET eVo|M, a, m n, sys) eos (2)= HeumG a.a, n: s. Bo Bc) 


Qu 7d* 5 -E+ all -ô) 
Hg. rm -ü42-0-E Bin 7 B*2-6-6 contrainte à, --n2d*2-0-eü-0*e-2-n 


Ko =]. Ôp m =2-0 &jg 72-6 


IT 


l-e 
8 Z m M " . | ] 
Pa = 0-7) f f z) Pa ()- Gun VD M, M 77h. B, mYr mrm ) er dur (z)- HeunG, la, nm Čr, m Ên mY nmn mZ | 
Qa, t =q+2-y-6-a(7+6-2) 
Qr nm =Q+3-y-ô-£=2-ß 


B, a =8+3-y-8-e=2-ã H z 
Ban - P E contrainte à, ,,,2-n-2- B 2-n B=n+2 
G=y+0+e-B-1=y+0+e-n-3 


2m -2-ü-n*5-y-ó-& 


Q? 


lina 2-y Ó, Iul -2-ó Ermm = 2-E 


l-e 
LI) y 1-6 Z DLH LIL HI m z n DLILIL( \ M n a 
En (z ) =Z ( SZ ) í z ] P € ) —Q € VoM, (a Ogg 77", Doa Yrm Bon ] et Ps (z ) - HeunG, [n v(i nm O ara Bru MED Ó ua z | 


Voici quelques exemples de polynómes de degré 1 : 


| acp [Boat +5) 


Haly -8)- B -1-2a— Jai B 


Ti 
z 
T 
= 
N 
— 
<< 
o, 
m 
D 


B-1-2a + Val 


> 
Il 
= 


> 
Il 
= 


1+a(2+7-6) -4ayl+ B z 


1c a(2 y -8)f -4ayl1+ B 2 
rist- ft a rit- (o at 


qo =q -ay(l-8) q =q -ay(1- 6) 


) £—y-2 Je alı B) y 3 4aylı B 2 
) £—y 24 6 al B) y 3) 4aylı B 2 


—£ l-e 
x (i x at) at) [12) (eat) 
a ay a ay 


qd, 745 -Yü-s) a -a5 -Y(-5) 


atf at 8-4)ey -B- ir - Beat 6-4) -4a(y -2hr «5 -2- B) 
o Jat oi eereo- der Beirat eo - 9] atr - 2 o-2- 3] 
Resa) feat) Res) i] 
dy ell y) a(y +8 2) di — dir ell y) a(y +8 2) 
ati! or e—4 a(y B) Ve Ee—4 a(y B) 4a(y »" &€-2 8) 
ín e peer Bede dr BJ] -sar hr s 72-8) 
seat) iun) rohi) (can) 
qg =q" -(-s)-(i-yXt*aó) aq =q" -(1-e)-(1-7X1+a6) 
at - e Bao - B)-A(B -e (8 -5)) -4«(8 -248 -:+2 o) 
at" = e B alë B) (6 E af 5) Aa(8 28 £42 2 


II IH 


Qo = dio 


1 
qi = a gta 


ay 


RM) 0-5) (oz 


-y(-e*a(1-8)) 4 = 


ILI ILI 2 1-ó X i X nu 
0 | Fi (x)=(1 x) í J [n aB- 2) | 


gi — y(1- e+a(1- à)) 


8-4)- (y * e-4 


a(y 


5-4) - 4a(y - 2v 


ó 


ó 


i 
gl = 1 gd 


aly 


ô 4)+ (y &£—4 


l-e 
pi Z 1-7 (1 1-5 1 =) 1 X 
(est 6-5) i 


LIII , 


a(y+5-2)+y7+e 


do — io 


2 di 


IBIBIII 


du 


aly 


LIT, I 


dio 


8-4) -4aly -2)(y 


| PE (xx Dal Jut X i 
a aa) 


+a(y+8-2)+7+e-2 


Notation équivalente des polynómes de Heun selon Ronveaux 


Voici un tableau des notations équivalentes que l'on peut retrouver dans la littérature et dans ce 
document : 


Classes | Classes Terme de F-|Terme de F-| Terme de F- a B y à £ 
Ronveaux | du Homotopie 1 | Homotopie | Homotopie 3 

document 2 

courant 
I 0 1 1 1 -n B y 8 £ 
II I x 1 1 1+a-y=-n 1+B-y 2-y à £ 
III Il 1 (1-x)'? 1 1+a-ô=-n 1-8 y 2-8 £ 
IV II 1 1 (1-x/a)^* 1--0-e--n 14 - € y ò 2-& 
V LII x" (1-x)'? 1 2+a-y-8=-n | 2-p-y-8 2-y 2-5 £ 
VI LII x 1 (1-x/a)'* 2-a-y-g--n | 24 p-y-e 2-y à 2-g 
VII ILI 1 (1-x)'? (1-x/a)'* 2+a-5-e=-n | 2--8-& y 2-5 2-€ 
VIII LILIIT x" (1-x)'? (1-x/a)'* 2-B=-n 2-a 2-y 2-5 2g 


Liens avec la construction des polynômes de Lamé 


L'équation de Lamé est un cas particulier de l'équation de Heun. Pour cela rappelons la forme 
Jacobienne de l'équation de Heun, en posant z = sn°(u,a) et a =1/0 : 


(«ay - etra) ps lodea) (as y: nan (a) = 0 


Avec les valeurs yzó-£71/2, il vient : 


Ga) eof sw (u.a) -4^q)y(u) - 0 | 
Posons 4aB =-v(v +1) 2= -4a?q > y'(u)+(a-v(v «Ya? sn^(u,o)y(u) - 0 


En posant maintenant le paramètre a = -n/2, et n-2N de manière à obtenir un polynôme de degré 
N, il vient : 


4a =-v(v+1 
ue l ) a eS DE An es 
2 2 d 
g--L--N 
2 


Il 
© 


> y"(u)+(2-n(n+ osi (us o) Qu) 


D'aprés la construction de la récurrence et des valeurs propres de la matrice, il vient pour le 
polynôme de Heun/Lamé de type I (ou 0) : 


1 n c n+l _ 2N+1 l EET. zai 
> 2 p aer a 3 A--4Aaq et x, =( nn 
r,=40°R,=(2j+2)2j+1) 

q,=4a°Q,=4j(1+a°) 

40° P, - 4a (j 1e 5-1+ B)» 4a? (5-1 sp; 14 


p; =4@°P, 2a? (j -1- NJQj -2*2N «1)- -a?(QN - 2j  2(2N +2; -1) 


x,=0 x«-1 


M 2a (j-1- NJOj-2*2N +1) Dr ME 


lg, n 0 * m 0 | 
Dd ñ 0 . 0 
0 0 T 0 
X S [oXy xn] M,-.. 0 pys dya y3 0 » (M, -AL,)X, =0 
0 Pyra dy5 ya O 
Ooan 0 Pra dya "ya 
0 Pu 4x] 


Aym valeurs propres de My )= RC B DE LN M Ay aE d € (0,,2..., N } 


Ce qui est bien la récurrence du type I des polynómes de Lamé, donnée dans l'ouvrage de 1962 des 
auteurs F.M.Arscott, I.M.Khabaza « Tables of Lamé Polynomials » (voir Table I, ligne (1), page xvii). 
et dans l'article de 1965 de J.H.Wilkinson « The calculation of Lame polynomials ». 


Le type II (ou I) de polynôme de Heun/Lamé, se construit ainsi : 


Posant À--4a'q et x,-(-lfe, n=2N+1 


if 
ace) zs n l ses b aai 2N 
a H N=a+->a 
i z | DEL 2 2 
H, cvs bebas contrainte «v 2n-22N41 4àf =-v(v+1)=-(2N +1J2N +2) 
3 l EE E cuu S 
5-7 Me 2 2 


r,=4@R,=(2j+2Y2j+3) 


2 2 170 =l 
q; -4a (e ge a] (o? ; o +1)+ s) +1) ( Fa ) a E US. 


p, - 4a? P, - 4a (j 1e à (5-14 B) 2 4a (j-1 Mn) a’ (2N-2j+2X2N+2j+1) 


Ay, valeurs propres de M x (a )J2 iz) EXE AL! m e {0,1,2,.., N} 


Ce qui est exactement la récurrence donnée dans l'ouvrage des auteurs F.M.Arscott, I.M.Khabaza 
« Tables of Lamé Polynomials » (voir Table I, ligne (2), page xvii). 


Le type IlI (ou Il) de polynôme de Heun/Lamé, se construit ainsi : 


Posant A--4a^q et x, =(-1)'c, n-2N-«l 


1 
US ARTE goce qal ipsc 
H, ,-üt- Bu B contrainte «v 2n-22N41 4àfi =-v(v+1)=-(2N +1)2N +2) 
"n. 17 173 2 
r, - Ad? R, =(2j+2)2;+1) 
2 
a posed prora 
Qt n EJ x0 xol 
FX ja t (s; E ah, +P,X ja = 0 


a 3 : 
q; -axo, +2) =(25+1 «4a? j? 


p; - 4a? P, - Aa? (j 1+8, j 1+ 8, )= a’ QN-2j«2Y2N «2j«1) 


jzN 
Av, valeurs propres de M,(a)— Pj, (z) 2 N1-z? Y dogs (A. yj m € (0,,2,.., N} 
j-0 


Ce qui est bien la récurrence donnée dans l'ouvrage des auteurs F.M.Arscott, I. M.Khabaza « Tables 
of Lamé Polynomials » (voir Table I, ligne (3), page xvii). 


Le type IV (ou Ill) de polynôme de Heun/Lamé, se construit ainsi : 


Posant A--4a^q et x, 2 (- Ife; n-2N«l 


1 
du, Qin N &+> > à zx 
el "E M y 
H, Gy ct, Bu -B* contrainte «v 2n22N41 4àf =-v(v+1)=-(2N &«1(2N +2) 
"m 2 3 
1 1 3 =N+1=8$,=N+- 
Yu 2 Ôr 72 £y x B Pu 2 
r, =4æ°R, =(2j+2\2j+1) 
2 
1 1 1 
em IDasclelee;ze ja 
ož jlleen n Pas] 
1 à T x, =0 x-1 
«(o i-i +40 j | 
d /^ 4 rx, tla; A px, =0 
; M ^ . ; 1 
P =(j-1+8nj-1+ 8n (1 DEL 
p, -4d? P, 2a? (QN -2j  2/9N «2j 1) 


J=N 
at, valeurs propresde Myla)= P! (z)- ar Y x yala mef12..N] 
j-0 


Là encore c'est bien la récurrence donnée dans l'ouvrage des auteurs F.M.Arscott, I.M.Khabaza 
« Tables of Lamé Polynomials » (voir Table I, ligne (4), page xvii]. 


Le type V (ou Lil) de polynôme de Heun/Lamé, se construit ainsi : 


Posant A--4am'q et x,-(-l/e, n=2N 


1 
Qj — REMIS 


in 7-N41od 412a =-N 


Hj, Gin = a+ Bin = B +1 contrainte 4v=n=2N 4a = v(v+1)= 2NQN +1) 
3 3 1 > l x 3 
Vin 2 re, En Pt cue 


r, z 4o? n, =(2j+2(2j+3 
1 


) 
= [je] 


2 2 
sa'o as a? " jua) jare œ+4j(j+2)  (x,-0 x,=1 
5X (a, = ax, +PX;1= 0 
q; -u lo ere |-Qj tla’ -a°+4j(j+2)+a° +4 


> q,=(2;+1) a +(25+2) 
p,=4aP=-(2N-2jf2N+2j+1) 


n,m 


4e valeurs propres de M, (o )2 p^ (z) -zNl-z Qus +. n) j m € (0,2, ,N] 


On retrouve la récurrence donnée dans l'ouvrage des auteurs F.M.Arscott, I.M.Khabaza « Tables of 
Lamé Polynomials » (voir Table I, ligne (6), page xvii). 


Le type VI (ou LIII) de polynôme de Heun/Lamé, se construit ainsi : 


Posant À--4a/q et x,- (- le, n-2N 


finir Gym =-N+1=&+1>&=-N 
Hs 4ã,m=@&+1  BPim=B+l contrainte 4v 2-n-2N 48B --v(v«1)- -2N(2N +1) 
3 1 3 z 1 2 3 
Vin = 5 Ô; i = 2 Eim — 2 B Ms ? Prim N« 2 


| = 
1 AM 
e, =al + ira’) K jaeren- sai 
rx +(g,-Akx.+p,x., =0 
— 4, -4a7Q, «4a? «1- Qj «1j +a°(2j+2) žm tla =h rin 
p, 4a? P, - Ad (j VG, Jj 1 B, )e ^ QN -2j ON 2] 1) 


EN 
Ar» valeurs propres de M, (a) P) = zyl- Q?Z’ XM ed yj m € {0,1,2,.., N} 
j= 


n,m 


On retrouve la récurrence donnée dans l'ouvrage des auteurs F.M.Arscott, I.M.Khabaza « Tables of 
Lamé Polynomials » (voir Table I, ligne (7), page xvii). 


Le type VII (ou II,III) de polynôme de Heun/Lamé, se construit ainsi : 
Posant À--4aq et x,- (- le; n-2N 


= ss) 
Qui = 4 4 nc 


-N+1=8+1>Q=-N 


Čr m 


Hs (Uu 4l  Bym=B+l contrainte dv=n=2N 4GB=-v(v+1)}=-2N(2N +1) 
1 3 3 = 1 y 3 
Kaga 2 Ou = 2 €i gp — 2 B=N+ 2 Bim N+ 2 


r, 4a? R, 2 Qj«2)2j +1) 
4a?Q, - Af «a? Kj 1)9 a, 24a? 0, 3e a? - (Là 2j 1 
Pj =4a°P, = 4a" (jM, s Xj -1+ Pom )=-a QN -2)9N «2 1) 


j-N 
ALU valeurs propres de M (a) P" (z)2 41-z? J-a?z? D EM (ans yj m € {0,1,2,.., N} 
j-0 


n,m 


1 =0 


On retrouve la récurrence donnée dans l'ouvrage des auteurs F.M.Arscott, I.M.Khabaza « Tables of 


Lamé Polynomials » (voir Table I, ligne (5), page xvii). 
Le type VIII (ou I,II,.IIl) de polynôme de Heun/Lamé, se construit ainsi : 


Posant 2=-4a°q et x, 2(-1'e, n=2N+1 


lrum 74t l+a 


Õim 2 Qnm =-N+1 2-8 > B N+1 
Biri zm yv-n-2N«l 4àf = v(v + 1)= (2N } IN } 2) 
Has \G=y+ô+e-B-1-y+0+6-n-3 CO"trainte og --(2N +1) 


> inm =n+5-y-ô-E€ 


3 
Ô, nm = 2 


" 1 ~ 2 5 
a= Nat) > Bin m 72-0 NM 


Vin = 


3 
2 Er Ip III E 

r,=4æ°R, 2 Qj 212j 3) 

4a?Q, - Aj(j 2o? +1)> q, - 4679, (1o? )- (1? J25 2) 

p, - 4o? P, - 4a (j 1*8, s Vj 1 B, )- - QN -25J0N 2j 3) 


EN 
valeurs propres de M, (a) po (z)» zdl- z? J1-a?z? DE pes kj me (0,,2,.., N} 
j-0 


x;,-0 x-1 
ru, +l, -i) tpx,,-0 


JT IH 


PEE 


n,m n,m 


On retrouve la récurrence donnée dans l'ouvrage des auteurs F.M.Arscott, I.M.Khabaza « Tables of 


Lamé Polynomials » (voir Table I, ligne (8), page xvii). 


Wronskien des solutions de l'équation de Heun 


Il est possible de calculer la forme générale du Wronskien de deux solutions indépendantes de 
l'équation de Heun avec un méme jeu de paramétre. Prenons l'équation de Heun sous la forme 
suivante : 


Considérons deux solutions indépendantes, alors nous obtenons l'équation différentielle du 
Wronskien de l'équation de Heun : 


e L e (z — lr (z m ay Wronskienly, (z), Y3 (2) 0 


Z 
Il vient donc : Wronskien{y, (z), »,(z}}= C(q,a,a,B,7,6) ARE (z-1)° (z-a)* 
Supposons maintenant que l'on connaít explicitement une solution de premiére espéce de 
l'équation de Heun. On peut alors déterminer formellement la forme intégrale de la solution de 


deuxiéme espéce par la formule bien connue du Wronskien : 


Wronskienly, (z), » G)} = C(q,a,a, B, y,6) zo (z — D (z — a)* 


y.(z) 

dsl = zov (2) v.n) E Wronskieniy,(z), v,(z)} 
dE x (z) y (z) 

ON CIEN n Wronskien{y (z), y,(z)] 
PON jac] aa | 


— yz)» xz) M 


Cas particulier : solution triviale de l'équation de Heun : 


D. cte pod 


LA] Y mes 
L'équation de Heun : |? e+ per 2cğ z-z-a) 0 se simplifie lorsque le terme 


Contrainte de Fuchs @ + B +1=7+0+e 
&Bz=Q s'annule soit lorsque üB-0 et q=0 L'équation différentielle devient alors une équation 
du premier degré moyennant : (z) U ô d | '(z) 2 0 
Z 


l-z a-z 


dLog(y'(z)) _ EEE TE €» Log(y'(z))=-y Log(z)- 9 Log(1-z)-e Log(a — z)- C, 


Z 
vient: — |y(z)- 6 fatt” (170) *(a-0)* «C. 
—1iGüB-0 et q-0 


G+B+1=y+0+e 


En posant a=1/K° il vient la forme suivante de la solution dítes « triviale » : 


y(z)- 6 [at TLC 


k € [0,1] 
GB=0 et q=0 


G+B+1=y+0+e 


Cas particulier : solution hypergéométrique de l'équation de Heun 


Ici je présente certains liens permettant de passer de l'équation de Heun : 


e + > + E Eo- 0 HeunG,(a,4,&, B. y, 5:2) 


à l'équation hypergéométrique : 


Ash) lez (a eet) G)-a ba) - 0er) EEEN QE (De > Fais) 


z-1) 


z-1 


Dans certains cas particulier l'équation de Heun devient un équation hypergéométrique. Par 
exemple lorsque q-a6a et £=0, il vient : 


er 2e À y(z)= um pou cni (y +6)}y'(z)+ ap y(z)=0 


&+B+1=7+6 


C'est bien une équation hypergéométrique et étant donné que la solution de Fróbenius développée 
autour de z=0 est identique à la solution hypergéométrique ;F;, à une constante près, on a : 


&+B+1=y+06 > HeunG,(a,& B a,ã, By. + B Y-yiz e F (a. Bsy:z) 


Zt dz = a dz 
~ a 
"(z)4. 4. € '(z ap z)= = x 7] Æ 
y ERES )+ La” ) ds y" (z)« Lezipeeu Gs 
AIEO PE ez(-z)v()-zi-B)G)-ag »(2)-0 
ä+B+1=y+e 
Soit : Hemala Pa. By sra fr) 
a 
On a également lorsque y=0 et q-0 : 
e, z-1 
zu dz = (a -1)dz 
3B a-l 
" E NEST a ne b L 4B 
vetere peut ans vepres a 0-0 
Contrainte de Fuchs à - B *1-2y*óc& ez-z»G.-z(-B«))»y(G)«ag y2)-0 
&+B+1=8+e 


Soit : HeunG, (a,0, Q, B,0,6: 225 


R 
Ww 
N 
l 
ee 


D'autres valeurs de paramètres conduisent à des relations triviales. Par exemple si a=1 et q-a6 : 


De méme pour a-0 et q-0 et alors 


e=0=6-&+p-7+1 


»'G) ! t La rye 1 0-0 HeunG,(a,q,&,B,7,8:z)= Fg. Byz) 


Développement de solutions de l'équation de Heun en fonctions hypergéométriques 


La premiére mention d'un développement possible est donnée par N.Svartholm en 1939 dans son 
article rédigé en allemand : « Die Lósung der Fuchsschen Differentialgleichung zweiter Ordnung 
durch hypergeometrische Polynome ». En fin d'article l'auteur propose un développement de 
solution à l'équation de Heun sous la forme : 


ylz)=} c, ,F(o vv +8 +y;y;z) 


Quelques années plus tard en 1944 A.Erdélyi publie dans son article « CERTAIN EXPANSIONS OF 
SOLUTIONS OF THE HEUN EQUATION » une étude plus systématique des développements de la 
forme: 


j=% 


ylz)= Xe, ,F(A* j, u- j;7;2) 


j=0 


Je vais donner les principaux résultats établis par A.Erdélyi et les relier au développement proposé 
par N.Svartholm. Les fonctions de base du développement d'Erdélyi sont au nombre de 6 
catégories : 


P()=(-1) du. FA ju- yz) 


DARE T) 
P)- T(A«j) T(I-y «A^ j) 
1 Tü-usj) T(y-u j) 
P'(z)-(1-z)" ,F(I-8* A j1-8* u- 2-012) 
ig PU IEU- p 6) 
T(I-A-u*2j) 
Ja TA-u22j) 
P n res ere 


27 ,F(lI-y A* j]-y* u- j2-y;z) 


FL u- iyl-z) 


rayr(- E +4); E * TT CP 
Rer (erg ^ tried 
T(i-u) (-u); 

T(1- uJr(y - u) t- u),r( - u), 


EVA y+À+ jl-À taie) 


ET . . "m wr : . «T 
ea jl-ytu-jl-À*u yl ea jl-ytu-jl-À*u 2 


On forme toute combinaison linéaire de ces 6 types de fonctions : P(z "X T, P!(z 


Ces six catégories de fonctions et leurs combinaisons linéaires quelconque respectent l'équation 
différentielle hypergéométrique suivante : 


Z 


ô ; ; 
«G-n [peo [23 eoa Ju- Dr. ()=0 
Remarque : toutes les fonctions P.(z) font partie des 24 solutions données par Kummer pour 
l'équation hypergéométrique. 


Pour cette méme fonction l'équation fonctionnelle suivante est vérifiée (ce qui est le plus difficile à 
prouver en terme de calcul) : 


a+B+l=y+ô+e 


no z)- Ka P. (2) (o; - a) P (2) M, P.) 


où les trois coefficients K,LM sont indépendants de l'argument z et seulement fonctions des 
paramètres a,6,y,6,g,A,u et j : 


pU tü-hu Ij*8-u-1fi*r-u-1Xi-u) 

2 2j+4-u-1{2j+1-u Y | | 

"rc («& +B + X +4), (j-&+a]5-B+1{5-7+1Xj-u) 
ner a) ^ ruens) 


"  (-&sAsMG- Brasil -y A12) 
au (2j+2-u+1(2j+2-u+2) 


On peut aussi écrire l'équation fonctionnelle sous une forme légèrement plus « condensée » : 


ez(-1)P(2)+@ B-(2+ )(u- 3) P) Ka Pa) ^, P) M, Pale) 
he u-Ij* B-u-ii*r-u-1G-u) 


: (2j*A-u-12j*-A-4u-2) 
a LUs&- ni B-ujiev-uXi*2) G-&*AX;-Be2|u-rAXi- n) 
g (2j+4-u2j+4-u+1) (2j+4-u\2j+2-u-1) 
ERU &+2+1Y; B+a+1{5-7+4+1X5+2) 
d (2j-A-u41(2j-A-u2) 


G+B+1=7+0+e 


Alors tout développement de cette forme vérifie une loi de récurrence à trois termes, ainsi que sa 
condition initiale et une contrainte supplémentaire, comme suit : 


Kc, HL; -a)e; +M c 


J 


0 


R 


Comme donnée plus haut les expressions littérales des trois coefficients K,L,M sont les suivantes : 


1 
(«à -uYj* B-ujiev-uYj*2)  G-&*2Xi- B Aj-v +a- n) 
(2j-A-uJ2j-A-u-1l) (2j-A-uX2j*-A-u-1) 

y Lcd ki Bano A2) 
j (2j+4-u+1\2j+4-u+2) 


Démonstration partielle pour l'établissement de la récurrence 


Nous allons supposer que la relation fonctionnelle entre la dérivée première de la fonction et la 
fonction P'(z)P(z) et ses plus proches voisins P (z)})P (z) 4 lieu. Du développement 
j=% 


ylz)= c, P.(z) appliquons l'équation différentielle des fonctions de hypergéométriques servant de 


ER 


base: | (z-1) ine. ô jeep eu nn) ainsi que l'équation fonctionnelle : 
e z(e- P) B z-a- (o DG DG- a) PEK a Pale) (L, a) P(2) M, Pal) 


dans l'équation différentielle de Heun : 


vet EE * X y) du »&)- à 


ede -Ie-a) C) (ee 04 3 E |y à z-)r)-o 


Ire eoe S «ea 0 - 070) 
+z(z-1)e y'(e)+@ B z-a-(e- aX Nu 7)»(2)-0 


s e-o- ro 3] ro] euro] 


© 


+z(z-1)e P(2)+(@ B z-a-(z-aYA jXu- j)) PG) 


Il vient : 


Il 
© 


eYee-ne rea B z-4-a+ lu Me-a)nG)-0 


oramg e z (2-1) P («lg B z-q-(4+ j)(u-j)(z-a)}P(2)=K „a P „(z)+ (L, - a) P ^M, | P (z ) 
es (K a P. (z)« (L, - a) (C)+ M, P,.(z)-0 


Terme P,(z )O K;c;, * (L, -q)c, +M,c 


Ff 


Comme. s- 0| re Le dene Xa - 0n) 


Cela implique donc la nullité des termes de chaque fonction de Base P; et donc la relation de 
récurrence sur les coefficients cj: Kc +(L, - q)e; +M je, 70- Cette relation de récurrence est 


complétée par la condition d'observation de l'équation de Heun: soit 
Xl (k n Paz )+(z, -q)P.(z )-M,, Pc )= 0= M ,P,(:)=0 et le départ de la récurrence en j-0 devient 
alors © Terme P(z) 5 (4, - q)e +M, c, > (D, - 4) c, * Mc, =0- 

Vous l'avez compris le plus difficile est de montrer qu'a lieu la relation fonctionnelle : 


e z (z-1)P'(2)+ @ B z-a- (A D(u- D)G-2) BG) - Ka PG) (b - a) PG) M s Pi) 


C'est ce que nous allons démontrer enfin de secteur concernant les développements d'Erdélyi et de 
Svartholm. 


La contrainte MP ,(z)=0 Se traduit par une condition sur les paramètres : Atysytdtleaxpos 


Dans certains cas c'est le coefficient M, qui s'annule dans d'autre c'est la fonction P ;(z) qui est 
identiquement nulle : 


P{e)=( n zh Luskyz) | Pi) Cy qeu y+4,2=y+m;2-7;2) 


pilej- GITEA rg Luekyl-) — PG)-(- 2" R3 12-9*p2-61-:) 


T(- a) Py - u-1 


T(A-1]r(-y +4) ii; 1 ; I(-u-2) 1 
P (z): = 7 F] A-1,- +4;4-u-1;- Pie) ——27"",#%| u+l,2-y+u3-A+ u; 
wee tm JE 


Cela conduit aux deux situations suivantes : 


Situation |: f zie 4-P-*. dans ce cas effectivement : do (à ay Bl2 2-1) _ o et 
A-B  u=à-e | (A-u-14 - u) 
toutes les fonctions Pm conviennent 


A=y+68-1 de 

Situation 11:}}=Y7  H=È-I 
=ô u=y-l1 
A=1 H-ytó-2 


Lorsque ; 2y 45-1 u=0 : 


Seules trois fonctions P., sont identiquement nulles, alors que le coefficient M; ne s'annule pas : 


Lorsque =y „=8-1 alors la condition y, - (G-api- Bp -v- | o est respectée et toutes 
- 
(4- u -1XA - u) 
les fonctions Pm conviennent de facto 


Lorsque 3-1  y-ysó-2 alors la condition y — (.-aa- BA -YYA-1) | o est respectée et 
(4-u-1 A - u) 


toutes les fonctions Pm conviennent également 


Lorsque 1-5 u-y-L Seules trois fonctions P, sont identiquement nulles, alors que le 
P(z)&0  Pjz)e0 Pi(2)#0 
Pi(z)=0  Pi)eo  Pi(z)e0 


coefficient M., ne s'annule pas : 


Lien avec les travaux de Svartholm 


On retrouve la récurrence introduite quelques années auparavant par Svartholm en 1939 dans la 


construction des séries de type Il. Pour cela placons nous dans le cas 3-,}5-1 ji 0: 


a 
I 


Svartholm expose en réalité une série de la forme:y .— (j)- & Æ(-jj+y+6-Ly;z) qUi 
j-0 


coincide (à une constante multiplicative prés) avec la série d'Erdélyi de la forme DNO ye, P'(z) 
rdélyi 


jja * +(1, -qk, tM, = 
Kg «(E -ay, «Mg, -o 


j j-l j T 


Dans ce cas les relations de récurrence peuvent alors s'écrire : lt 


e. & jT(y *ó* j-1 nas 
Nissan) = Yoani l2) = c. ,EC j,j+y+ô-l;y;z )= c,( 17 (+ E ) g( Lj*y*6-hy;z) 
Comme j=0 p T(1- j) 
is T(y 4-14 j) 
.2(-1/ 
vut I(l+j) ^ 


Il vient donc : 


Tfy+6+j-2) à T(y+6+j-1) > T(y+8+j) = 
K L M 1) K, 1) L: 1} M. =0 
jja th; jC; * jÉj T ( ) r(j) E Î ( ) T(1+ j) €; j ( ) (2^ j) Cia J 
S r(j)  r(y«ó«j-1). y 1307-2 
d ^T(y--64j-2) T(+j) 3 j 
> L, =L, 
E T(2+j) r(y«8*j-1).- esce 
/ IT(y+6+j)  r(j) ^ytó«j-l 
Ce qui conduit à la récurrence suivante : 
g Usa BA ey M9 5-2) 
(2j+y+ô Jon y3 t 3) 
P ceni y, Usa Burke, (i-e y 6 - Bey o-i eo 1); 
j 7-8 JN 3 JH : | 
(2j+2X2j+2+1) (2j - AY2.j- A -1) 
g s cgo Bey opis ote 5) 


4 (2j+y+82j+y+8+1) 


Avec la substitution d'indice j->v et le changement global de signe opéré par l'auteur et q-h selon 
les notations de Svartholm, il vient : 


g (rero very e6-2Yv-à-1fv - B-1) 


i Qv -y -6-2)2v y 6-3) 
(v - av B kv y Xv - A) (v-àcy6-Mv- Beyeo-1Jve6-0)7 


Vaid QE ME ver ren Qv « AYov «A -1) 

rese - a eye SÀv- Byso) 

i (2v+y+82v+y+8+1) 

Nv+8-1)=(-1Xv+7+6-2)+à B 

ee Y d 
(v -a)v* Bkv-yXyv«2) (v-à^ 

Comme 4av(v+y+8-1) ADIT 

[Gti eroi ehe rie ern g za (4728-02 À) 

g rrr 8-2 ey e5-2)«8 B) 

: Qv*y48-3Qv «y48-2) 

(= t se Dogs DEP), à aversa A] 


(v«à- 
(v-&+7y 
+515 +7+6-1{v+6-1); - 
Qv - A2v  A-1) 


EU. Qv cy *8-22v «y +8) 


wav re yo) À) 
j Qv +7 «82v +7 +6 +1) 


K c, SHE, e he, e, -0 


J 


Convergence des séries d'Erdélyi 


Je vais essayer maintenant de brosser rapidement les critéres de la convergence de ces séries de 
fonctions hypergéométriques, selon ce que je peux en lire dans l'article d'Erdélyi : « CERTAIN 
EXPANSIONS OF SOLUTIONS OF THE HEUN EQUATION » et ce que je peux en tester numériquement 
avec Mathematica. 


Développement de type Il: | 


> 
II 
T Q 


Les développements avec les fonctions hypergéométrique P'(z) P'(z) P'(), P*(z) P*(z) convergent tous à 


deux conditions : Premiére condition : q est la racine de l'équation transcendantale de la fraction 
continue : 


K,/M, » 
L,[M;- 


-IL/M,- 


on peut aussi dire que q est la racine de l'équation issue de l'annulation d'un determinant infini 
suivant : 


LM O0 zm Q^ uu. 
K L M 0 D ume 
0  … T 0 t 0 
M = 0 K, L,, M, OÙ ss 
j4 j-l JA 
M jet (, F ak, tK, 0 K L M 0 (M-g1)C=0 
J J J 
0 Kj, Li, Mj 
C= ese nec 


<q tel que det(M = gl) =0 


Deuxième condition: l'argument z dans le plan complexe est à l'intérieur de l'ellipse donc les 
uo sont z=0 et z=1 et qui passe par le point z=a. Cette dernière A a = équation : 


|’ ul . La condition pour l'argument z s'écrit dienen: Ó Z Is 
mm 1+ l1- l4 Ri 14 
Z a 


Le développement avec p'(;) n'est jamais convergent quelque soit la valeur de q partout dans le 
J 


plan complexe 


Le développement avec P*(;)-est convergent quelque soit la valeur de q pour les valeurs de z en 
J 


Le iyi 
dehors de l'ellipse, soit zl a, Typiquement pour z-x»a sur l'axe des réels, le 


Be fuel 
Z a 


développement pi p*(z) est convergent. 
J X 


j-0 


A-y 
Développement de type ll : A-cy | gu-ó-l 

A-ó 

VES 


Pour 2-7  u=8-1, les développements Se, P) Se pt) sont convergents sur z€[0,1] 
dud ? JT 
j-0 j-0 


lorsque q est la racine de l'équation transcendantale précédente. 


; j=% j=% je 
Pour 3-5  u=7-1, les développements Ye, P(z) Ye, Pc) Ye, P‘(;) sont convergents sur 
j-0 j-0 j-0 
z€[0,1] lorsque q est la racine de l'équation transcendantale précédente. La convergence pour le 
développement Se, p*(z)a lieu car la fonction ,« p*(z) est un polynóme de degré j en z. 
J^. J 
j-0 


Pour 1-1 u=y7+6-2, les développements Ns P?(z) Se, p*(z) sont convergents sur z€[0,1] 
Jo d AU 
j-0 j-0 


lorsque q est la racine de l'équation transcendantale précédente. 


2 je j=% [ee 
Pour 3=5+y-1 | u0, les développements Sc, P(z) Ze P*(z), Ze p*(z) Sont convergents 


j=0 


sur z€[0,1] lorsque q est la racine de l'équation transcendantale précédente. La convergence pour 


le développement y c, P‘(z) a lieu car la fonction z“ p*(z) est un polynóme de degré j en z. 
D LEAN; J 


j-0 


| 1 | 1 
1-./1 1-41 
En théorie tous ces cas de convergence ont lieu à l'intérieur de l'ellipse Z> a, mais 


l+ pr 1+ jd 
a 


Z 


numériquement je ne suis pas parvenu à le vérifier, aussi je laisse l'intervalle z€[0,1] pour plus de 

sûreté. En dehors de ces cas de figure il n'y a pas de convergence en dehors de l'ellipse pour aucune 

des fonctions P! (z), P? (z), P’ (z) P*(z) P*(z) P$(z) que la valeur q soit contrainte à une racine de 
J cat Tog tr T Tu 


l'équation transcendantale ou non. 


Quelques graphes de développements hypergéométriques, cas de convergence 


Commentaire : lorsque la valeur propre q est définie comme l'une des racines de l'équation 
transcendantale, soit par annulation du déterminant, ou par détermination des valeurs propres 
d'un systéme linéaire, soit par formation des fractions continues associées et résolution de 
l'équation transcendantale par l'algorithme de Bouwkamp par exemple, alors les coefficients du 
développement sont déterminés par les vecteurs propres de la matrice associée. Dans ce cas la 
normalisation co-1 n'est pas forcément respectée mais la solution approchée obtenue l'est à une 
constante multiplicative prés. 


Dans la mesure du possible j'ai donc vérifié de maniére empirique la convergence des 
développements proposés par A.Erdélyi, dans l'article original de 1944 « CERTAIN EXPANSIONS OF 
SOLUTIONS OF THE HEUN EQUATION ». Les cas de convergence y sont explicités en pages 66 à 68. 


Convergence de la série Se, P(z.surzf[O1i]avec ; 25 ^... 
F9 


Paramètres : a=1.25 a= 1/3 62 1/3 y2 1/3 ó-2/3£za*6*41- y- 6 ÀA- a u-6-&€ 
q--0.0250046 est la première valeur propre de la matrice tridiagonale M : 


—— HeunErdelyiExpansion(pos, a, a, B, y, Ô, À, p, Z, 1) 
—— HeunErdelyiExpansionPrime(pos, a, a, B, y, 6, À, p, z, 1) 
—— HeunErdelyiExpansionPrime2(pos, a, a, B, y, 6, À, p, Z, 1) 


Convergence de la série Se, P (z)-sur z€[0,1] avec } -g usH-s 
j-0 


Paramètres : a=1.25 a = 1/3 62 1/3 y 1/3 622/3 £za*6*1- y- ő Â= au=86-E€ 
q--0.0250046 est la première valeur propre de la matrice tridiagonale M : 


—— HeunErdelyiExpansion(pos, a, a, B, y, 6, À, p, Z, 2) 
—— HeunErdelyiExpansionPrime(pos, a, a, B, y, 6, À, p, Z, 2) 
—— HeunErdelyiExpansionPrime2(pos, a, a, B, y, 6, À, p, Z, 2) 


Convergence de la série So PXz;)surzf[O1i]avec 32g  p=8-e 
j=0 


Paramètres : a=1.25 a = 1/3 6= 1/3 y= 1/3 622/3£za*6*1- y- ő À= au=86-E€ 
q=-0.0250046 est la première valeur propre de la matrice tridiagonale M : 


—— HeunErdelyiExpansion(pos, a, æ, B, y, 6, À, H, Z, 3) 
—— HeunErdelyiExpansionPrime(pos, a, a, B, y, 6, À, p, Z, 3) 
—— HeunErdelyiExpansionPrime2(pos, a, a, B, y, 6, À, p, Z, 3) 


Convergence de la série jy P:(-) sur z€]0,1] avec ; .; u=ğ-e 
j-0 


Paramètres : a=1.25 a = 2/7 6= 7/13 y2 2/3 = 1/5 £za* 64 1- y- ó À- a u-6-&€ 
q-0.0472702 est la premiére valeur propre de la matrice tridiagonale M : 


—— HeunkErdelyiExpansion(pos, a, a, B, y, 6, À, p, Z, 4) 


—— HeunErdelyiExpansionPrime(pos, a, a, B, y, 6, À, H, Z, 4) 
—— HeunErdelyiExpansionPrime2(pos, a, a, B, y, 6, À, p, Z, 4) 


;, . dem z 
Convergence de la série y, p*(z) sur z£la,+ee] avec ; y u=ğ-e 
j=0 
Paramètres : a=1.25 a = 2/7 6= 7/13 y=2/3 = 1/5 £za* 64 1- y- 6 À=q@u=68- e 


q=-3.21918 est la troisième valeur propre de la matrice tridiagonale M : 


0.006 


0.004 F 


t€ —— HeunErdelyiExpansion(pos, a, a, B, y, 6, À, p, z, 5) 
—— HeunErdelyiExpansionPrime(pos, a, a, B, y, 6, À, p, z, 5) 
- 0.002 —— HeunErdelyiExpansionPrime2(pos, a, a, B, y, 6, À, p, z, 5) 


-0.004 


- 0.006 H 


-0.008 t 


p. de 
Convergence de la série Sc, p*(z).sur z£[a,+ce] avec » -g D m 
F 4J 
j-0 


Paramètres : a-1.2a = 2/3 62 3/7 y2 1/3 ó22/7 £za* 6 1- y- ó Â= au=6-E€# 
quelque soit la valeur de q, ici q=3.23, avec 40 termes pour le développement, résolution directe de 
la récurrence donc co=1 : 


N 


—— HeunErdelyiExpansion(nTerm, a, q, a, B, y, 6, À, p, z, 5) 
—— HeunErdelyiExpansionPrime(nTerm, a, q, a, B, y, 6, À, u, z, 5) 


—— HeunErdelyiExpansionPrime2(nTerm, a, q, a, B, y, 6, À, p, z, 5) 


NE 


z dX 
Convergence de la série Sc, P'(z)-SUr z£]O,1] avec ;-, 28-1 
j-0 


Paramètres :a=2  a-7/8 4623/4 y=3/11 86=2/7 e=a+6+1-y-ô À=y u=6-1 
g=-15.7153 est la quatrième valeur propre de la matrice tridiagonale M : 


— HeunErdelyiExpansion(pos, a, a, B, y, 6, À, pi, z, TypeFun) 


o —— HeunErdelyiExpansionPrime(pos, a, æ, B, y, 6, À, p, z, TypeFun) 
—— HeunErdelyiExpansionPrime2(pos, a, a, B, y, 6, À, p, z, TypeFun) 


Convergence de la série Se P*(£) sur z£]0,1] avec ;-, | ,-28-1 
j-0 


Paramétres:a-1.2  a-27/8 6=3/4 y-1/3 8=2/3 e=a+6+1-y-ô ÀA-2y u=6-1 
q--2.706668 est la troisième valeur propre de la matrice tridiagonale M : 


6r 


N 


—— HeunErdelyiExpansion(pos, a, a, B, y, 6, À, p, z, TypeFun) 
—— HeunErdelyiExpansionPrime(pos, a, a, B, y, 6, À, p, z, TypeFun) 


—— HeunErdelyiExpansionPrime2(pos, a, a, B, y, 6, À, p, z, TypeFun) 


Convergence de la série STe P? (z) -sur z€[0,1] avec ;-| u=y+68-2 
j=0 


Paramètres : a=1.2 &=7/8 6=3/4 y-1/3 &=2/3 e=a+6+1-y-ô À=y u=6-1 
q=-1.92554 est la deuxième valeur propre de la matrice tridiagonale M : 


—— HeunErdelyiExpansion(pos, a, a, B, y, 6, À, p, z, TypeFun) 
—— HeunErdelyiExpansionPrime(pos, a, a, B, y, 6, À, p, z, TypeFun) 


—— HeunErdelyiExpansionPrime2(pos, a, a, B, y, 6, À, p, z, TypeFun) 


Convergence de la série Se, P*(z)-sur z€[0,1] avec ;-| u=y+68-2 
j=0 


Paramètres : a=1.2 &=7/8 6=3/4 y=3/5 6212/13 £=qa+ 6+ 1-y-6 A-y u=6-1 
q=0.0379612 est la première valeur propre de la matrice tridiagonale M : 


—— HeunErdelyiExpansion(pos, a, a, B, y, 6, À, p, z, TypeFun) 
—— HeunErdelyiExpansionPrime(pos, a, a, B, y, 6, À, p, z, TypeFun) 
—— HeunErdelyiExpansionPrime2(pos, a, a, B, y, 6, À, p, z, TypeFun) 


Convergence de la série p P?(-) sur z€l0,1]avec 325 ^ ,-y-1 


j-0 
Paramètres :a=2  a-7/8 4623/4 y=1/3 6=2/9 e=a+6+1-y-ô Azó g-y-1 
q-0.255647 est la première valeur propre de la matrice tridiagonale M : 


3 


— HeunErdelyiExpansion(pos, a, a, B, y, 6, À, u, z, TypeFun) 
—— HeunErdelyiExpansionPrime(pos, a, a, B, y, 6, À, p, z, TypeFun) 
—— HeunErdelyiExpansionPrime2(pos, a, a, B, y, 6, À, p, z, TypeFun) 


Convergence de la série 37 P'(-) sur z€]0,1] avec 325 u=y-1 
j-0 


Paramètres :a=2  a-7/8 4623/4 y=1/3 6=2/9 e=a+6+1-y-ô Azó u=y-1 
q--8.15284 est la troisième valeur propre de la matrice tridiagonale M : 


—— HeunErdelyiExpansion(pos, a, a, B, y, 6, À, p, z, TypeFun) 


—— HeunErdelyiExpansionPrime(pos, a, a, B, y, 6, À, p, z, TypeFun) 
—— HeunErdelyiExpansionPrime2(pos, a, a, B, y, 6, À, p, z, TypeFun) 


Convergence de la série Se, P'(z)-sur z€[0,1] avec ;-,.5-1 u=0 
j=0 


Paramètres :a=2  a-7/8 6=3/4 y=7/13 6=2/27 £€=qæ+ 64 1-y-ó À=y+ő-1 u=0 
q=0.578384 est la première valeur propre de la matrice tridiagonale M : 


—— HeunErdelyiExpansion(pos, a, æ, B, y, 6, À, p, z, TypeFun) 
—— HeunErdelyiExpansionPrime(pos, a, a, B, y, 6, À, p, z, TypeFun) 
—— HeunErdelyiExpansionPrime2(pos, a, a, B, y, 6, À, u, z, TypeFun) 


Convergence de la série Sa P*(z).sur ze[o1]avec ;-,.5-1 u=0 
j-0 


Paramètres :a=2  a-7/8 6=3/4 y=7/13 6=2/27 e=a+6+1-y-ô A-y4ó-1 u=0 
q--0.4616 est la deuxième valeur propre de la matrice tridiagonale M : 


—— HeunErdelyiExpansion(pos, a, a, B, y, 6, À, p, z, TypeFun) 
—— HeunErdelyiExpansionPrime(pos, a, a, B, y, 6, À, p, z, TypeFun) 
—— HeunErdelyiExpansionPrime2(pos, a, a, B, y, 6, À, p, z, TypeFun) 


-10r 


Convergence de la série ye, p*(z).sur ze[o1]avec ;-y.5-1  yu-0 
Jd 


j-0 
Dans ce cas la fonction hypergéométrique dans la définition de Ps est un polynôme de degré 
maximal j en 1/z : 


I(A-u*2j) iai | : | | 
Pé(z)- z"B5 AÀl-y*u-jil-A-u-2j; 
1) ere een Ca Run ee - 
A-ytó-l ; „1 , u 
jd HN -jl-y-j2-y-ó e Polynóme de degré j 


Donc la fonction eni ND gf jl-y-j2-y-ó ui Polynóme de degré j est elle- 
T(L+ (y + 7) z 


méme un polynóme de degré maximal j. C'est l'une des conditions de convergence de la série 


j=% , s z . 

Sc, FE (z) avec celle que a soit racine de l'équation transcendantale. 
d d 

j=0 


Paramètres : a=2 q= 1/3 6= 1/3 y= 1/3 = 2/7 £& za4641-y-ó6 À= y+ő-1 u=0 
q=-0.829572 est la deuxième valeur propre de la matrice tridiagonale M : 


— HeunErdelyiExpansion(pos, a, a, B, y, 6, À, p, z, TypeFun) 
— HeunErdelyiExpansionPrime(pos, a, a, B, y, 6, À, u, z, TypeFun) 


-10r —— HeunErdelyiExpansionPrime2(pos, a, a, B, y, 6, À, u, z, TypeFun) 


Convergence de la série Se, PS(£)-sur z£]0,1] avec ;-5 u=y-1 


j-0 
Dans ce cas la fonction hypergéométrique dans la définition de Ps est un polynôme de degré 
maximal j en 1/z : 


T(A-u2j) i | | 2 
P*(z)- uj p 1 dur 
(z) Tü-u«jr(-usj). ,QÉZ|H-JÀA-Ytu-]p H ^7 


ÀA-ó 1 
| — rr l- j,-j;,y 0 —-2j; | Polynóme de degré j 
u-y-l z 


Donc la fonction ;«p(;)— T6 «1-7 2j) ZGU NE. 2 Polynómededegrà: j estelle: 
TQ -y « ri 5) z 


méme un polynóme de degré maximal j. C'est l'une des conditions de convergence de la série 


[ze ; 5 z P 

Ye, P‘(;) avec celle que a soit racine de l'équation transcendantale. 
J J 

j=0 


Paramètres :a=2  a-1/3 6=7/9 y=3/11 6=2/7 e=a+6+1-y-6ô À=6 yu-y-1 
q=-0.164481 est la première valeur propre de la matrice tridiagonale M : 


—— HeunErdelyiExpansion(pos, a, a, B, y, 6, À, p, z, TypeFun) 
—— HeunErdelyiExpansionPrime(pos, a, a, B, y, 6, À, u, z, TypeFun) 
—— HeunErdelyiExpansionPrime2(pos, a, a, B, y, 6, À, p, z, TypeFun) 


Remarque : Il est clair que les développements d'Erdelyi sont également possibles avec toutes les 
formes de Kummer, puisque chacune de ces formes sont des combinaisons linéaires des deux 
premiéres fonctions d'Erdélyi. Et que moyennant un bon choix de constante multiplicative cette 
combinaison linéaire sera indépendante de l'indice j. Dans ce cas on peut légitimement se poser la 
question d'un développement possible avec une fonction hyper-géométrique s'annulant en 0 et en 
1, mais hélas les formes de Kummer ne présente pas de telles fonctions. Il suffit de regarder les 
formules (15.10.12) et (15.10.14) données dans « NIST Handbook of Mathematical Functions », qui 
sont les seuls forment régulières en z=0 et z-1. Les produits de la forme z^*(1-z)'*^ dans w2(z) ou 
ou w4(z) correspondent à des formes équivalentes de ces fonctions qui ne s'annulent chacune 
qu'en un seul point z-0 ou z-1. 


Quelques relations fonctionnelles concernant les fonctions hypergéométriques 


Prenons la première des fonctions de base P; utilisée par Erdélyi toujours à une constante 
multiplicative prés : 


P. (EER (A+ j+u-j-1;7;2) 
P (HE+ j-Lu-j+l;y;2) 


On peut formaliser le problème en établissant deux relations aux « plus proches voisins » de la 
forme: 


PEE ju- j;yz)> 


se- GR bsez]- Aae) „F(a-1,b+1;c;z)+B(a,b,c) Æla+1b-1;c;z)+ Cla,b,c) F(a,b;c;z) 


Z F (a,b; c;2)= Dla,b,c) „F(a-1,b+1;c;z)+ E(a,b,c) „Fla+l,b-1;c;z)+E(a,b,c) F(a,b;c;z) 


sle- lebscz]--d Æla-1b;c;z)+(c-a-bz) Fabi) 


(b-a) Ela,b;c;z)+a,F(a+1,b;c;z)-b,F (a,b+1;c;z)=0 
5(b-a«1) F(a-15;cz)4a-1) Ela,b;c;z)-b,Fla-1,b+1;c;z)=0 


De plus : 
b-a+ 
2) Fabia) ch Fette ET : F(a;cz) 
JC WT NE hot , js m )740,0, C; 


En tenant compte de la symétrie a«-»b des fonctions hyper-géométriques : 


d -c)b 
Z E-licsFlabiaz)- a 


De méme : choses) = 
z -4+ 


Fett E: Fla,b;cz) 


Alba-heshi( E 


: z Æba;c;z) 


MEAE -e)a lets) : : Fab; c;2) 
dz +] +] 


z ina Eus on 2+), Flab:c:2) 
+] b-a+l a-b4 


-—— i b) 
PUE L — 7 — aebl-1) Flabie 
ab d b-a+l "E oz b- 2 IJa- nj“ + e), (a, 61) 


Soit finalement la première relation : 

ü-cC F( 

— 3 Fa-bbesez P = Frthb-1;62) 
la-bfb-a+1) ^ (b- A bl c 

La seconde relation s'établit à l'aide de l'expression suivante : 


Pinen = F(a,b-1;c;z) 
a 


E ot E RUE 


xr ee US 
dz mE zi 


noire j 


SuF(ab cz Fab;cz) E 


He-Micz)-— Æla-1,b-1; 2) 


Or (b-a+1),E(a-1,b;c;z)+a-1},F (a,b ;c;z)-b,F(a-1,b+1;c;z)=0 


b a-l 
Abel baies) lab:62) 


b-l 
EE o — à mr Fn: 
De plus ,F(a-15- laa) e- Frs = F(ab-1;c2) 
D'où 2,F (ab: ;Gz EFI lab;c;z He e e-ble- JE eh pisos) 
a- Ifb-a) a b-a 
2uF(ab cz EF, lie. F(a- Wie Æ(a,b-1; 2) 
a 
c-d 


exar Rieb;ez) e-Micz) Alab-lica 


eb? 


Ce qui donne la deuxiéme relation recherchée et établie ainsi ce que l'on souhaite : 


Pob; cz E. F(nb;c) 
dP(a,b;c;z) 2c-a-b-] ü-c 
-1 IVO rnt " b P b; PES 
i2 c b-a+ifa-b+1) i Dre cem 
abeu v er) b(c-a) 
a-b+1)b-a+1) (b-a+1fa-b) 


Pla-1b+1;c2)+ P(a+1,b-1;c;z) 


b-c 
(b-aa-b+1) 


P(a-1,b+1;c;z)+ aet) ») Pla+1b-l;c2) 


l- boa] 


z P(a,b; c;2)= 


La deuxième fonction de base P; est obtenue par la transformation de Kummer suivante : 


LFob;c;z) 927. Fa- 
Transformation de Kummer a-c«1291-y*4* j 
b-c+l>l-y+u-j 


cdb-cd;2-61) Ti- ; 
-J*^t] 


Tea) 


hu 


„Ell-y+4+jl-y+4-j;2-7;2) 


Prenons alors la première et la deuxième relation en l'appliquant à la fonction 


Fla-es1b-e41;2-c;) il vient : 


is vostre 


dz 


zz" Fla-c+lb-c+1;2-c;2)=4+ 


De plus ale- JEE Rlar hbc? 
z 


Il vient finalement : 


-e;z)}= 


a-l 
la-bjb-a«1) 
b-l 
la-bla-b+1) ° 
a+b-| 
p-a fab 


a-c+1)(b-c+1)- 


2 Fla-c+1b-c+1;2-c:2}+ 


z™ Fla-c+2,b-c;2-c2)- 


Æla-c,b-c+2;,2-c2)- 
Fla-c+2,b-c;2-c2)- 


Fla-c+1b-c+1;2-c;7) 


2 Fla-cb-c+2;2-c;2) 


-(1-c)z * E(a-c*15-e41;2-c;:) 


Lab cla+b 1) a +b Ia c*lfb c«l) 


t- 


(a-b+Ib-a+1)  (b-a+ifa-b+1) 


zl JEP c+l,b-c+1;2 ezz} 


l 
{eh 


U-ba-b+1) 


(a-b+1fb-a 


Ann 


(b-c 
t-il 


la-5)b-a«1) 


aber) i. e c+1b-c+1:2 c;z) bloet) ve Pla e43 b-e:2-e;zht 
(a-b+1fb-a+1) la-b+1]b-0) 


(-c)z2,Fla-c+1b-c+#1:2-c;2)- 


Piel) ef Rl-er-etli2-cià] 
Z 


( ds cHb-cH;2-c;z) 
) t Ih csi ezl) = Fla-c+2,b-c:2-c:2)+ 
| | Fa-erhh-ei2-eis) 


4yt- fruits 


all-afb-c+1) 


b-a+1]b-0) 


2 Fla-ch-c+2:2-c1) 


Ce qui donne donc les deux relations pour la fonction de base : ‘Pla, bic: DE 2 Fla-c+1b-c+1;2-c;2) : 
b(i-b\a-c+1) 
la-b+1)b-0) a 

HD pa :c;z) ped core Lb+l;c:2) 


b-a+15-0) ° 


Pla+1 nesl- ai Lbsl;c;z) 


(b-a+1)b-0) 
| 


Pla b:c2) 


dPabic;z)_ abLc- la+b+1) 
n God) 


P 


is 
METRE 


z Pla,b;c;z)= 


2ab- 2ab-c(a+b-1) 1) 
(a-b+1fb-a+] 


La troisième fonction de base P; obtenue par une transformation de Kummer s'écrit : 


Hakeai) y mhr 


| FU ji- j;81- 
A+j+u-j+l-y=A+u+l-y=a+b+l-y=6 | RS ar LA tk jio 


Transformation de Kummer | 


Il s'agit donc d'une fonction de la forme L Fla b:cl- z} pour lequel a lieu les deux relations 
375.923 


suivantes : 


d od b-c poU a+b+1-2c aa 
ck- pr F(a,b; cC- scii ET x Ela -1,b+1;c;¢) (a-bla-b+1) Fab Le) G-0+1%0-5+) Han] 


Ce em ent Fab T an VE Rest) EL Æl-1btl;cé) 


(-l-z et Plab;c;2k,Rlabic;l-2)= 
dP{a,b;c:2) c-(a+b+1) a-c 
l-z] — =ab b; —— — Pla-l,b+1;c;z)- 
p zs za eben ge c eene 
c(a«b-1-à -9 +1 
a-b+1fb-a+t) 


b-c 
——— Pla+l,b-1;c; 
(a-bja-b41) 4 f ) 


olimga: «et Pla+1b-l;c;z)+ M J Pla- Lb+1; 651) 


(a-b+1fb-a) ` 


Pour la fonction de base P,, par le changement z->1-z, on déduit les relations suivante : 


z P(a,b;c;z)= 


Pla,b;c;z)=(1-z}",Ffa-c+1,b-c¢+1;2-¢;1-7) 

re dlabic;z)_abla+b+1}-26) TEN b(l-bfa-c+1) bla cl), y dl-0)b-c+l) ME 
eJ M are re Le er aa Mets 
M NELLE) , cen) T "a cHJ(i- d), des 
Hehe C risas] ere er LL et jer 


Pour les fonctions P; et Ps, il est plus aisé de calculer directement la relation fonctionnelle avec les 
paramétres données par Erdélyi. 


Relations fonctionnelles d'Erdelyi sur les fonctions hyper-géométriques de base 


A partir des calculs qui précédent il s'agit d'établir maintenant la relation fonctionnelle qu'Erdelyi 
introduit dans son article de 1944, et de former la combinaison linéaire suivante entre les deux 
expressions p'(z)z P(z): ez(z- P/'()«(« B- (A j)(u-) z PG) Kj Pj) ^, PG) M, P, G)- 


Jj ju 


Fonction P; : sachant que : 


P(z)«C, ,F(A* ju- j:y:z)- C1 


T(4+ j) "HN ĉi 
= A+ ju- jiy) >E , === 
r(i-u+ j)? ( | C, A+j Ca j-u 


et que les relations aux plus proches voisins ont été établies : 


P(a,b;c;z}H,F (a,b;c;z) 
x dlabicz) àc-a-b-] 7 RE " b-c A 
dn Aa LEM ur arre 1 Len rm 
pen) do) eu a 
Hebe aai PP n tene eng tne 


s b-u-j c=y>a-b=4-u+2j b-a=u-4-2j a4bzÀA*uzy40-lza*f-c 
BLAN" ] p), Ue i-rhsi- D p (yl lui LNH) p k 
do (u-A-2j+1f4-u+2j+1) Fe re ul  (u-4-2jfA-u+2j+1) ^ 
À VE EU 
AE TC er 77 PR CT TC ^ TPS PE ^ 
ep P(e)=K PACA PMP) e=a+8-4-u 
(i-u+a+Blu-j)-a8_(i+a-u)(+8-ufi+r-uli+1-u) 
(i-u+2j+u-à-2j) Qj«eA-uQjÀ-u) 
Lebe iu DOG Du i) D) o a ia BA nor Au) 
^ (4-u+2j+1\u-4-2j+1) 
Ljsa- nk B-nVier- uli à) U-a sii Pre afin) 
(2j+4-u\2j+4-u+) (2j+4-u\2j+4-u-1) 
Jw (i-re) Auc sa DAUBA- +A- 
aem uem a NR E E 


K={7-u+jfl+j-u) 


Fonction P?; : sachant que : 


l-y+4+ j) , | 6; Et ®© __ née 
AGE C, z ” Ell- -ytA* jl-y*u-j;2-y; z)- (yore i) F(I-y 44* jI-y*u-j;2-y;z)9 E - —= - 
Hr-u*j) '' C, l-y+4+j Cj l-y+4-j 


et que les relations aux plus proches voisins ont été établies : 


Plab;c;z)=2",Fla-ctib-c+l;2-c;2) 
jeter) abe lebe] c bei) alizal- H) i- a aci! PNE 
d la-b+1)b-a+1) ^ b-a+1)b-a) ^ la-b+1fa-b) 
l-a 


o) aL; a UD) pit 17 
(b-a+1fa-b) a-b 


2ab-c(a+b-1) (ub; je cH) 


Gba) 77 


z Pla,b;c;z)= 


Il vient : 


a=4+j b=u-j czy2a-bzA-u42j b-a=u-1-2j a*bzAtyuzasf-e 

Cie D CoA] 

Ci l-y+4+j Ca l-y*u-j 

MP) lu rn) ac i Ire (y, Un ar) p 

z- - | D pf;}+ P? (z)+ PA) 
d  (-utj#4fu-1-2j#) ! "tj- IsA-uoj-u«4) © (Aà-a2j a -n2j) * 

P Ir Ira) so, ore a, (Catia) s 
(A-us2j«lu-A-2j41) ^^ Qj-1sA-nfA- 2j) "^ -u2j Mf -u22j) ! 


| dP'(z) 


Ba PE PA PEN B 


scs Pis) ran P- uli r- Xi 1-n) 
ORALE T TP (2j+1+4-u\2j+4-u) 

cvs asa n eR) -e s PAX re Aint 
Macr dpa bj- 
Seu Dr oet) 8-0 Du DOG la )- 7 T) 
i (u-A-2j+1fA-u+2j+1) 
sani B-nVje- uli 3) (j-a+4hj-B+4Nj-y+4j-u) 

Qjsi- aii ue) Qjei- uj ca -u-1) 


Fonction P2; : Sachant que 


SUCRE E | YN; 
= J 
(l-u+ jy -u+ j) (5 j-1 Xj -y*4) 


— 


T(i« jjrü-y«4- j Sra C, 
l SE ju- j;ô;l-z7)> LEE - - 
r(i-u+jE(y-u+ j)? Ca (+jl-y+4+j) Ca G-uly-u+j-)) 


P(z-C, .K(A* j,u- j;5;1-z)- 


et que les relations aux plus proches voisins ont été établies : 


Plab;c:2k,Rlabic;1-2)= 
PA LAZULON ab 2c- (a1) C 


(caius 


Plab;c;z)+ 


P{a-1,b+1:c:2)+ ————— 
p b-aon uet i 


alc-b) 


(a-b+1}a-b) 


la+b-1)-a°-b +1 


Pla+lb-l:c; 
(a-b+1fb-a+1) Ua 


ZEE d nk 
pap nie 


zPlab;cz)=< Plab;c;z}+ 


Il vient : 

a=À+j b=u-j c=0=a-b=1-u+2j b-azu-À-2j atb=l+u=y+0-1l=a+p-e 

A+u=7+0-1506=1+u+1-7 

DEAN dci a*f-e«1-20 PRE Wire) p tft) p | 
dz (u-A-2j+1f4-u+2j+1) ^^ (A-u+2jfA-u+2j-1)  (1-u+2jfA-u+2j+) ' 

pg. onore Hift p, BA pp erh n 


(a-u+2j+1{u-4-2j+1) S (u-4-2j4u-2-2)) ^ A 
dP'{z 
ea) apte me 2 eh bar e 


osos à 08-Uxitellu-ÿ) G+a-n)i+8-uli+7-uli+1-u) 
OR AE) Aa THIVEP DNI 
pen 288-579) i-a - Pea -r Aiit 
Masa rd ee Davos) 
p ue p-ei728)a B i) Deer) 8-61) 0j] if 
j (u- 4-2 A - ui2j 1) 
Uta ki B- kie -uli3) U-a i PAX afin) 
Qj«i-ufjsà-ul) Qj«A-ufojeA-u-1 


Fonction P^; : sachant que : 


P'()-(-:)* (o j41-8,4- j 1-8;2-5;1 jo har 


et que les relations aux plus proches voisins ont été établies : 
Pla b:e:z)=(1-z}",Fla-e+l,b-¢+1;2-¢;1-2) 
(jetez) aber} 18) pi piepel- T ! cte) 
ko (a-bHfb-ad) b-a+1fa-b) (a-b+1fb-a) 


r Plae; D pipe Bd ee ug e) 
UO gebe CC 777 e ^ ebela | 


Pla+1b-l;c;2) 


Il vient : 


ü-À*j b=u-j czó2a-b-A-u42j b-a=u-À-2j atb=1+u=y+0-1=a+p-e 
Atu=7+0-156=1+u+1-y7 A+l-ô=y-u us-ózy-À 
cay Go e pet-n8] » PRAE) CN) p Uc P) 
& (U-uHj4fu-2-2j4l) ^'^ (2j-1+4-ufh-u+2j) "^ Q-aau-a ? 
M M ; f ; ; 
p) ej ui) p) W+1-8- ci) p) nr) aj) P'(; 
annaba bete te) 


eset apte me, ea bar 


osos à 88-9 isgu-i) _(j+a-u)lj+b-ulj+y-uj+-u) 
DT TE Qaia 
Lu ns a 8B elei) G-oAg- Pri rA 4-1) 
Mati ba Gel 
Lea p-e+1-28)+la B-l i-o 9-141 jt] 
(u-2-2j+]4-u+2j+1) 
j*a- nj B-uXjr- n2) U-ae- P Aj Afin) 
(Qj+4-u\2j+4-u+) Qj«A- uj À-u-l) 


Passons aux fonctions P; et Ps, liées aux transformations de Kummer : 


Flab LE jo Ph 


l 
ucc 
Transformation de Kummer : 
„Elab LE j pu 


beta) 
Z 


Soit avec les paramètres donnés par Erdélyi : 


a=Àti b=u-j ezy 


| | 
PÍaj;c;z)- 2: ‘Alu c+l;a-b+]; j Pun -y+4+j;l+4-u+2j; : 
E 
4 n | | l 
Plabic;z)=z Ælbb-c+l;b-atli-|=2 Fl u-jl-y+u-j;l+u-4-2j;- 
1 Z 


Précisément ce sont les fonctions suivantes : 


T(i* jrü-y 4j) à | 2 
P(z)- | ElA+j1-y+4+;j:1+4 2j; 
"(z) core t Hd ML 


lü-u-2j) — 4 | | | J 
P*(z)- FL F| u-jl- -ji-A*u-2j;- 
“(z -aiya i UTJA- Y +H]; *u D7 


Pour trouver les relations fonctionnelles sur ces fonctions, le plus simple est encore d'utiliser les 
formules de connexion de Kummer concernant les diverses fonctions indépendantes et solutions de 
l'équation hyper-géométriques. D'après « NIST Handbook of Mathematical Functions Chapitre 15 
Hypergeometric Function », a lieu les combinaisons linéaires suivantes, ou encore formules de 
connexion : 


w(z)- Elab; c;z},E(a,b;c;z) 
m(z)=P(a,b;c;z)=z,Fla-c+1,b-¢+1;2-¢;2) 


v(e)-e"adsz)- snae 
Z 


codes duit J 


Dead) C, eue dt), ( 
(2) H Doy p) E Kc) 

c-lrb-a- 
wet? HOME a) w(z) 


Soit avec les paramétres et constantes multiplicatives données par Erdelyi : 
G-A*j b=u-j  czy 


wei uhr) Cy eat py 


w(z)- 7 FA j-y+Lu-j-y+1;2-y;z)=(-1} iC P; 


I(I-y 444 j) / 

— ili) FIA+iA+i-y+1:14+2j- sz EDS 0417 vous P UNE 
wj) € Z 2 ( tJ, t] y* , + J Ht Jt (A+ j)r-y 14 j) 7 
Eur eiua 

Tü-242j 


leen T(i-4-u42j) Pi) 
l ) 


-dVi 


Meter la ju-j-y+l;u-1-2j+1;- 


> Pee” th-y)Ple)+etry-)el) 


s dep). EL 2 Cr ET gea To a2 Eu — 
7 Toy-aejoprtu- j-r +i- jr j) ' Tü-u jr(u- ipi -2- jp 2 j) ' 
ARR a ART 
lcu ir- je TET =) zi Sus u) 


Or Neue MT 


T(y-u+ j|r(u j sn 1} TEE 


í 
ri-i-ü-reiej)-- T " 1-j) a U 1) 


us DSi (y = inl 4) ue ere Tür Sin u)Sinla (7-2) à: 
z Sin( (i- u) Pee x Sin(z (i- u) ^e) 


Pise 


Donc les fonctions P^; et P5; sont des combinaisons linéaires des fonctions P'; et P^; indépendantes 
de j. Et à ce titre les fonctions P”; et P5, vérifient automatiquement la relation fonctionnelle 
d'Erdélyi. C'est d'ailleurs bien ce qu'indique Erdélyi dans son article. J'aurais d'ailleurs pu également 
utiliser cette remarque fait par l'auteur pour démontrer la relation fonctionnelle en ce qui concerne 
les fonctions P?; et P^; . On a donc bien démontré la relation : 


e z(z-1P'()*la B-(4+ J)(u- j)z P(z)e Kja Pa) ^, P) M, Pa) 
x _G+a-u -A*B-u-M*r-u-YG-u) 
2 GT E CIE) 
A LG ra-uXi* B-uXivr-uXj*2) G-acAXji- BeAXi-y «AX n) 
: Qj-A-u)aj*A-u*1) Qj-A-u)2j*A-u-1) 
i (j-a*-A«1(j- B A«1(j- y X A4X1Xj - 2) 
f (2j+4-u+1X2j+2-u+2) 


a+B+1l=7+0+e 


A titre indicatif je donne les combinaisons linéaires associées aux fonctions P’; et P*; en fonctions 
de P et P^: 
w(z)- P(a,b We z)-; F (a,b cz) 
w(z)=P{abic;z)=2"*,Rla-c+1b-c+1;2-c:2) 
w(z)- (abc; z Fab a b-c1;1-z) 
)= 
)= 


Z 


Plab:c;z)}=(1-2)"",Flc-a,c-b;c-a-b+1;1-2) 
r(I-c ener), y, Pe Thiet eil), 
T(a-c«ir(b- c1) | T(aJr(b) , 


uU T(1-cJr(c- aba Tie-Ine-a-bt) o 


r(- a] (1-b) T(c-a)r(c- 5) 


G-A*tj b=u-j c=y  Atuzytó-l c-a-b+1=1+7-1-u=2-6 


wz 


T 


c-a=y-À-j=u-j-0+1 c-b=y-u+j=4+j-ô+1 a+b-c+l=ô 


nakabara Un 


T(y- uj) pi) 


=z" F(A ] D? ly 
w(z)-z (+) y*lu-j-y* y;z)= ( lr” 


T T'ü-u-* jriy-u*j 
wc EF ju j:0:-2)- vd ne 


wz)-(L-z)^ ,F(A* j41-6,u- j41-8;2-8;1-z)- P. jd 
vet e er p 
P()-Gy T -7) 8-8) P()« C1) i iro- » pe 


T-i- jas j) ^ T(y-A- jjü-y 4j) ? 


» OP sisi) 


F(-u- ju j)- C1) m 
Pi(z)+ 


p- U ln 6-4), 


p)- r(1- yr - 6 )Sin(r À) pg) r(y a in (y-4) p 


T T 


Développements finis en fonctions hypergéométriques d'Ishkhanian 


L'auteur A.M.Ishkhanian et ses condisciples ont proposé des développements finis en fonctions 
hypergéométriques des solutions de l'équation de Heun de paramétres : 


ee bre 59 geo 


z z-l z-a (z-1Xz-a 


dans une série d'articles: en 2005, A.M.Ishkhanyan, R.Sokhoyan, D.Melikdzanian « New 
hypergeometric series solutions to the general Heun equation », en 2014, T.A.lshkhanyan, 
A.M.Ishkhanyan « Hypergeometric expansions of the solutions of the general Heun equation ». Ces 
développements sont de deux types : 


Type (1) LOS Cj „Fla, Biy, — iz) 


Type (1) »(z)- Cj „F(a, B;y, + j:2) 


j=0 


Ces deux développements conduisent aux relations de récurrence à trois termes suivantes sur les 
coefficients c; : 


Type (I) Ale, (B -a)e, Cc, =0 Type (x) A ej (B -g)e, e cilc, =0 
PORC Yo -J 1)(B Vo jy Vo j+) 
A! = (a Wy, irs ntj 1) J yo *j-l 
BO zaa f-(a-1ky,-j-Yy c-r *j)*a(y-n*iXa*B-0n*j) KBP =aa B-(a-My e j-My *e-n - D)*a(- 5 - Ja B9 -J) 
co afe n*i*(-n*i*Ur-n*i) C -(a-1ky, e jy +E- -j-1) 
yo-j-l 


Les conditions d'arrét du développement en j-0 sont déterminées par l'annulation du terme C de la 
récurrence en j=-1 et conduisent à contraindre le paramètre y, à des valeurs précises, soit : 


Yo = 
Type (I) C0 -02(a-)(8-v(v-5)-09 170 = 8 
Vo =Y 
=] 
Type (II) ct -o- Ur +e-n)=02 
Vo EYE 


Dans les développements de type (Il), la valeur yo-1 n'est pas licite puisque la récurrence présente 
, ; Pun Ei 0 " " 1 
pour le terme c., une forme indéterminée : Type (11) 4 x La relation de récurrence n'est 


donc plus respectée. 


Pour que les développements soit finis, il convient que le terme A de la récurrence s'annule à son 
tour pour une valeur donnée N+1 de l'indice j : 


Type (I ) Atc 


FA 


«(81 -q)e, + Che, =0 Type (I) Ale 


Jal 


yy-a-»e-a-y-N e=a-y-N 
A =0=y+e-n+N=0= 17 =B—e=p-y-N  |A=0=y=y+e le=B-y-N 
Yo =y >E=-N E€=-N 


Il s'avère qu'en réalité les développements de type I et Il correspondent à des développements 
identiques. Comparons les deux développements finis de type II : 


Type (1) AUR. + (a'^ -qJe, " ce -0 


Type (11) 


AU, E 


j=N 


Yo = à s=a-y-N= y(z)- 


j “ji 


DN a; 


=y+e e=a-y-N= ylz) 


i" 


j- 


ak, teea =0 


For 


Les fonctions hypergéométriques sont donc les mêmes. 


Fla, B;a- j; F3 & Ra, Bia f -Niz)- T 2 N-j 
j-0 
F 
Fla, B; B- jz)- Y y F(o.B:B j -Niz)- T 2 Nj 
j-0 
jcN | v às 
Biy- Az- e Fo Biy*j-Niz)e- T 2 N-j 


j= =0 


SS UR (a, B;a  j - N;z) 


j-0 


(«Ye Rl. fpe j- 2) 


Pour les coefficients de la récurrence des développements de type I et Il: 


70€ 
e=a-y-N 


(n. 
C= 


P =(a-Ija-jXj-N-1) _ 
Baa --ie-j-M- tales) TE 
n _a(i+)(B-a+j#1)(7-a+ j*1) Vo =7 
dá a-j-l 
AP -(a-148— jj -N-1) n 
B” =aa f-(a-AKB- j-1(j- N)ea(v- B^ j Va j) erp 
o ala- Bj «UG B i1) h-r 
l B-j-1 

A” -(a-1ky - ij- N-1) 40 = 

--N 
(a- 1\y- ja Wi- N)+a j(a+B-y+j) i ETE pin 
a(a-y*j*(B-y* j*)(i) Wy -rtezy c- 
pgs 


() 
B J 
i y B'-a ab- 
£ = 


qo aN-j+)(B-a+N-j+I)(y-a+N-j+1) 
A a-N+j-1 
MET: B* j(a-1fa- N* j-1)-*a(y-a* N- jf N- j) 


i. o-ife-N« +) 


C S5 


40 afa B N- jN- jer - Be N- jt) 
: B-N«j-1 


BU =aa B+j(a-WB-N+j-1)+a(y-B+N-jfa+N-j) 
M z Xa-1K8- N jXj1) 


e 


ala-y+N-j+I)(8-y+N-j+I)(y-y+N-j+1) 
y-N*j-l 


-aa f j(a-My- N^ j-1)*a(N-jfa* B-y* N- j) 
-a-Wr- Niki) 


Pour chacune des trois situations de développement fini, il vient systématiquement les relations 


suivantes : 


(nz 
Aj^- 


cp, P=) 


0). 4). 
Cj = AN- 


Dans ces conditions les coefficients de la 


récurrence de type I ordonnés en sens inverse, suivent les mêmes relations de récurrence que la 
récurrence de type ll : 


Type (I) AP, + (8! 3 q), +CPE; 4 =0 
dc, He eat 


=> AN En 


Cy. ja 7 Cza 7 Oja 
uy (gt Ja (n 7: 
> Aj Cya tB; -4)6, ; tC; €. = 04, 76; 7€; 
Cy 


A une constante multiplicative prés, on peut donc écrire : 


Type (1) Type (11) n-2B-2»)-; & Flo 8:84 -N:z)- | e „Fla, B; B j- N;z) 


Dans ces conditions et conformément aux articles des auteurs | A.M.Ishkhanyan, R.Sokhoyan, 
D.Melikdzanian, T.A.Ishkhanyan je n'étudierais que les développements finis de type II. Précisons la 


forme des deux premiers développements de type Il : 


y _a(N-j#1)(B-a+N-j+1)(y-a+N-j+1) 
e=a-7-N Le um 
li E »e)= Yr, „F(a, b;a+ j- N;z) BU = a a B+ j (a IM N+j 1)+a (y a+N iXB+N j) 
' Cf? =-(a-1a -N+ jj +1) 
4n 4l -B+N-j+)(N-j+)(y-p+N-j+1) 
e=ß-y-N e debe 
li Bep y)- Me, F(a Bie j- Ni) 4B saa Be j(a-IB-N j-1)ea(y-B*N-jfa*N-j) 
^ c -(a-I(8-N e jj +1) 


Les deux premiers développements sont symétriques par l'inter-change des paramètres a<->6. Il 


suffit de regarder donc le premier développement : Live Æ(a,fia+j-N:z): Appliquons une 


j-0 

transformation connue des fonctions hypergéométriques * Fla bas jz)=(1-2),F(ja-b+ j;a+ j;z) 

fonctions de ce développement : 
j=N 


X22 Ye; (1-2) ^ RON + j.a- B-N* jja- N j;z) 


j-0 


T 
Or s+y+8=a+p+l=a+ô N31-8--N-f yz)- Ye, (i- zy" , EC N* j,a- B- N ja - N- j;z) 


j= 


[m j-N -p-N«ji l IEN-j i-N -p-N-*i l 

De plus EE E ja- TA nir A. ` un ue 

1-0 1 . I-0 l 0 
a ETC (a- B-N- j) z ylz) 

HN V e. (1-2 JT} I E 
Mis l ? Zel z) 1=0 (a-N+;j) L (p 


Polynome de degré < N 


Comme les articles indiquent que la fonction yG-z est elle-méme une solution d'une autre 


équation de Heun avec des paramètres différents. Dès lors les deux premiers développements ne 
sont jamais qu'une re-écriture des solutions polynomiales de Heun à l'aide de fonctions 
hypergéométriques. Ces deux développements ne forment donc pas de nouvelles solutions 


algébriques. 


Voici quand méme les expressions suffisamment simples et polynomiales des deux premiers 
développements finis (N=0 et N-1). 


Pour N-O, il vient : 


y= E=aQ-Y ó-1«f 2 qayf -ay(B-1) 


Pour N=1, il vient : 


N=1 20a £za-y-l ózacfsl-y-e-2 422 B 26-2 

BU =aa B+a(y-a+1{B+1)=a a B+ae(1-6) 

c” =(1-aa-1) 

q-B) — q-a(a B+eli-6) 

qu ^ fafa- 

q qa -1-a(f +y+26y))+apy(-a aL BL y) -0 

eq *qa-1-a(6 -2«y(28 -3)))- ay(6 -2-a +a(5-1{1+7))=0 


tar [EIE i-e [Esta 


(a -1) a- 


8st en rre SN) 


a-l (L-aa-1) 


mz "m + Cle, =0>c=c, 


ER 


q racine de 


Par contre le troisième développement des solutions de type Il est plus intéressant comme l'indique 
T.A.Ishkhanyan dans l'article de 2018 (je le note développement 11.3), car il ne conduit pas à une 
forme polynomiale, mais il reste fini : 


40 afr N- jy N- j81(N-j*1) 
i y-N*j-l 


jzN 
y&)- Xe, „Ela, By * j- N;z) BU =aag H jla ly N+)j 1)+ a(N jla tB-y+N j 
P Cz (a-1fy - N^ jKj-1) 


ged Ale -r+N-j+)(B-y+N-j+I)(N-j+1) 
i y-N*j-l 
Je le note donc : 159 -aa 8. j(a-1fy -N j-1)&a(N - ja * B- y N- 5j) €t les coefficients c; sont solutions 


J 


c") z -(a- 1 y - N jj) 


d'un système linéaire fini d'ordre N+1 : 


Ae «(gir? “dk, «Cle, =( 

RG D Go) Le es 0 ][«] Tol 
AU» gedag cu D aoo 0 a | |o 
D amb. dise rdi s. A 
un s ‘6 AD ms CES [lea] [0 

| 0 n e 9 0 AP a a L0] 


Dés lors que ce développement est fini, cela signifie que le déterminant associé du systéme linéaire 
des coefficients c; doit s'annuler, soit : 


[pe ceo 0 ut P i ME 
A pu 3 q cu 3) 0 " LS 0 
0 A9 B"? -q cu) 0 " 0 
Det m =0 
. o 4 og qu 
| 0 0 AD sg 


Ce qui détermine donc q comme racines d'un polynôme de degré N+1. Par exemple le cas le plus 
simple étant N-0 pour une seule équation linéaire, il vient trivialement : 


j202 BU? 249 g M =" -q] De(M)-0« q- B? = à a f 


Les premières expressions de ces développements finis sont assez triviales à déterminer. Pour N=0, 

il vient comme on vient de le voir: i e se P . | s'agit donc du cas de 
n=1 bec .R(o.fi.z) 

dégénérescence de l'équation de Heun vers l'équation hypergéométrique. 


Pour N=1, il vient : 


q racines de g°-q{i-y+ala+B+2a B)}+a a B(L-a 1» B)-y)- 0 
j=l 
Yo=Y-1 yz)- Ye; „Fla, By * j-5z)2 e, „Fla, B;y-1;z)+c, „Fla, yz) 
j=0 
B” - gk, +c =0 
uius i i. q-aa B*a(1-6) 
: CU  ? (l-afy-1) 


£--] 


BU -ag B-a(a B*1-y)-aa B «a(8-1) 


NP I IR DEC RS RS 
yz)-. Fla. B; l; )+ l-a -1) „F(a, fy; ) 


Equation du second degré 
q!-q(I-y*a(a* B«2a B)}+a a B(Le o1» B)- y)-0 
© (q-aa Ba(1-5) qaa B (a-11- y))- al - a Lea - 7 XL B - y) -0 


Pour N-2, il vient : 


q racines de 


E (a-aa B) +(q-a a B)(4a a(3*a* f)-24 y)*2a (a lj Be 


x(q aa p 2al1+ a + B) 2+2y)+2a(a 1q aa Pl+a1+6)=0 


j=2 
y) Sc, Fla, Bir+ji-kz)=c Fla, sy - 2:2) e, ,R(a,B:y-Lz)+e, ,F(a,Br:7) 
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yy 2y-2 


(Bt Ie dk, +e =0 


S T pestifera 

C z(1I- aky-2 1 © 09» 7 f[(rafy-2 

BU = a a B+2a (5-1) > E 2«q)-a(a B y +qla+B+2a B « 26 2j 

go _2la-1+2)(8-7+2) o ae Emm Pao p) 206-1) 26 -1fo-2) 
de à Aa-1} (y My -2) 


C? - X1 - aXy -1) 
B” -aa B«(a-1fy -2)«a(5-2) 
x(z-. Fa. iy -2;2)4 Na Fla, B:r-152)+ 
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«a (e Bl(1+a)1+ 8)+2(5 -1))+2(6 -1 8 - 2)) 
Xa -1] (y -1ky -2) 


a, b;y;z) 


Question de la convergence ou de la divergence des séries infinis de type (I) et (II) 


On a donc vu que les séries finies d'Ishkhanyan de type I sont en réalité des séries finies de type II. 
De plus les séries finies de type II.1 et II.2 sont apparentées aux polynómes de Heun. Il reste donc à 
se poser la question de la convergence des séries infinies de type 1.1, 1.2, 1.3 et II.3. 


Je ne donne pas de résultats formels en la matiére, mais j'ai testé numériquement les convergences 
des séries initiées par Ishkhanyan. La convergence des séries infinies de type 1.1, 1.2, 1.3 est vérifiée 
lorsque q est solution de l'équation déterminantale « infinie » suivante (ici je ne donne que le 
développement I.1 mais il en est de méme avec les développements 1.2 et 1.3) 


Lee. gu Oo n. ds ow 0 
At) gl) -q co» 0 
0 A) BU» -q ce» 0 . 
Lim Det E B m w s P " -0 


No 
ʻi » s 0 A Bex C4 
0 a i T 0 AO BU) q] 


En prenant une valeur de N suffisamment grande les premières racines q peuvent suffisamment 
bien approximées. A titre d'exemples prenons : a=1.2, a=7/8, 6=3/4, y=7/13 et 6=2/27. Sachant 
que £e=a+6+1-y-6 pour le cas yo-a, les graphes de la fonction et de ses deux premières dérivées 
sont : 


—— HeunlshkhanianExpansion(pos, a, a, B, y, Ó, €, y0, z, 1) 
—— HeunlshkhanianExpansionPrime(pos, a, a, B, y, 6, €, y0, z, 1) 


—— HeunlshkhanianExpansionPrime2(pos, a, a, B, y, Ó, €, y0, z, 1) 


0.1 | 0.2 0.3 04 0.5 


Prenons maintenant : a=1.2, a=7/8, 8=3/4, y=7/13 et 622/27. Sachant que £=a+8+1-y-6 pour le cas 
Vo76, les graphes de la fonction et de ses deux premières dérivées sont : 


—— HeunlshkhanianExpansion(pos, a, a, B, y, 6, €, y0, z, 1) 
60 —— HeunlshkhanianExpansionPrime(pos, a, a, B, y, Ó, €, y0, z, 1) 


—— HeunlshkhanianExpansionPrime2(pos, a, a, B, y, õ, €, y0, z, 1) 


Le graphe reprenant l'équation différentielle de Heun est « numériquement » nul : 


-1.5 x 10711 | 


Et finalement prenons : a=1.2, a-7/8, 6=3/4, y=7/13 et 6-2/27. Sachant que £=a+68+1-y-6 pour le 
cas yo-y, les graphes de la fonction et de ses deux premières dérivées sont : 


—— HeunlshkhanianExpansion(pos, a, a, B, y, Ó, €, y0, z, 1) 
—— HeunlshkhanianExpansionPrime(pos, a, a, B, y, Ó, €, y0, z, 1) 


—— HeunlshkhanianExpansionPrime2(pos, a, a, B, y, 6, €, y0, z, 1) 


Le graphe reprenant l'équation différentielle de Heun est « numériquement » nulle : 


Important : je n'ai pas trouvé de convergence numérique pour les séries infinies de type II (soit pour 
yo-€*y) méme lorsque q est racine de l'équation déterminantale « infinie »s suivante: 


pu) 2 cu) 0 
AUS E ce? 0 . 
0 A B? -q cm 0 
Lim Det 
N>% 
o 4 Wie C 
| 0 0 A3) pu) -q| 


Ces solutions étudiées de développement 1l.3 par les auteurs A.M.Ishkhanyan, R.Sokhoyan, 
D.Melikdzanian, T.A.lshkhanyan ont cette particularité que ce sont en fait des fonctions 
hypergéométriques généralisées. Elles ont été successivement exhibées dans divers articles : 


- 1994, J.Letessier, G.Valent J.Wimp « Some differential equations satisfied by hypergeometric 
functions » 

- 2011-2012, K.Takemura « Heun's equation, generalized hypergeometric function and exceptional 
Jacobi polynomial » 

- 2018, T.A.Ishkhanyan, A.M.Ishkhanyan, « Gauss-hypergeometric expansions of the general Heun 
functions governed by two-term recurrence relations » 

- 2018, T.A.lshkhanyan, T.A.Shahverdyan, A.M.Ishkhanyan, « Expansions of the Solutions of the 
General Heun Equation Governed by Two-Term Recurrence » 


Dans ces trois articles les notations et les expressions ne sont pas forcément toujours les mémes, 
aussi vais-je m'efforcais d'unifier les propos, autant que faire se peut. L'idée de base est de montrer 
que la construction de ces solutions particuliéres des fonctions de Heun à l'aide de fonctions 


hypergéométriques 2F1, conduit à des solutions déterminées par une récurrence à deux termes et 
sont en réalité elles-mémes des fonctions généralisées hypergéométriques. 


Commençons simplement par le rappel de l'équation différentielle hypergéométrique de la fonction 
2F1: 


»}F (a. f:r:2) 
Eee ee] (enl) stretto) ao 


H dz dz\ dz 


On a vu ci-avant que l'équation de Heun se réduisait à une telle équation dans le cas N=0. 


Les fonctions généralisées hypergéométriques se notent : (;). F. k x 
bb, b. 


13-7299 


Elles sont solutions de l'équation différentielle suivantes : f É La) E ri É d, b, ipa 
dz d 


iz 


Une des solutions de cette équation différentielle est la fonction généralisée hypergéométrique 
dont le développement en série autour du point singulier régulier z=0 est le suivant : 


»Xz,F, Lee ue 2» (a), OMS T(a, -1).r(a, +1) r(b).1(b,) zl 


bb] ES) A) f cm rar) r6 e.r(o, 1) n 


Dans l'article de K.Takemura, il note l'opérateur différentiel de l'équation des fonctions généralisées 
hypergéométrique : ; comme étant l'opérateur différentiel dont le coefficient de l'ordre 


4,0) s... p ibi sba sb, 
maximal est 1, soit l'opérateur différentiel « monomial » pour lequel la fonction généralisée 
hypergéométrique est solution de l'équation :, ,,,_ ,[y(z))-0 


Par exemple pour une fonction 3F2 (voir K.Takemura), il vient l'opérateur différentiel suivant: 


d,,d,,d 
Xe" : 5 


bb 
d rm ta 3) zb +b, +1) d | 
E dz z (z-1) d? 
Eai | (aa, +aa, +a,a, +a, +a,+a,+1)z-bb, d NE 
i z'(z -1) dz z'(z -1) 


En 1994, les auteurs J.Letessier, G.Valent, J.Wimp étudie l'équation différentielle des fonctions 
généralisées hypergéométriques de la forme suivante : 


4,, 0» ,..., d, € +1,..,e, +1 
4 
, 


EE b,,b,...b,e,..e, 


Ils montrent que ces fonctions satisfont à une équation différentielle dont les coefficients sont des 
polynômes. 


L'un des principaux algorithmes, exposé dans l'article « Some differential equations satisfied by 
hypergeometric functions », consiste à calculer l'équation différentielle d'une fonction dérivée 
donnée à partir de l'expression de l'équation différentielle de la fonction de base. Plus précisément 
soit l'équation différentielle de degré s de la fonction f : 


l=s 


nd ideo we (T Yat 


FSI 
1=0 dx 


coefficients p(x) sont calculés par l'algorithme suivant : 


d i: d! d Fs di | 
(2 -Yu9. CDD X(x)=ax+b c;=c+a j 


j=0 


CS RE i 


e()e ue) Ge)  je0.s+ oou) oue) n) 


first MONCET skal] 


k=l k=l 


Les coefficients px). p, (x) ont les expressions simplifiés suivantes : 


Lt) Eu (ed Ys DAL E H 5, o ()- es) 


E s-l 
e x) st ) Eat} SIE: DAC "m 
l k=s-1 5 ( j s-l k=s k=s 
o pul x)+ ut x DAMES J- Hst mal} Eat R(x) AOC DAS (0) 
k=0 z 


R(x)- Du DA) R'(x) z * (u, (x, (x)+ H (u'l) L'(x)= EL MOE a M 
Y Gh Rue 
ES X X ER. XIX JC, = EN X xh. v XJ AX eR) a v XJ AX 
ah G)-- x P x6 [ssh best eaeh -- 5093 eut) 
y E REY n 
Dep Futaba) e SRG RE st ER) 
cr x ARX ipd xJlAx)- Rix) R'ix)- u,'lx y xyl Ax 
saifa te TS Ro rs) Et] 
E! aR(x) $ V i ; ka 
= ne rose B e tok ett 
de) {Ru (St) 


a, Be. +1 


L'application de cet algorithme est par exemple applicable à la fonction ez, e| ; } puisque 


X8 


nous savons des fonctions généralisées hypergéométrique que : 


f E «proe "a 9.4 (z) 


à e 


Plus généralement si l'on construit la fonction g(z) de la manière suivante : 


„B.e +l, e, + 
doen E 2 ZA 


9.6.8, e, 


a, p,e *l..,e, +1 
fef | «| 
19356,5563 
Elle suit parfaitement le schéma algorithmique proposée dans l'article « Some differential 
equations satisfied by hypergeometric functions », sachant que l'équation différentielle est 
constamment du second degré. 


En l'occurrence l'algorithme a la forme suivante : 


k=l k=l 


po) RG)us Gr) S()u x) ntm 


PAT e cu a RD) 


Comme exemple ils donnent en fin d'article (expressions 3.22 et 3.23) la fonction : 


dE ea 


/,& 
qui sous certaines contraintes entre les paramètres g 5, , satisfait à l'équation de Heun. 


Nous pouvons donc calculer assez facilement l'expression de l'équation Mp c de la fonction 


g(x) à partir de celle de l'équation différentielle de la fonction f(x), soit f(x), s" a} en suivant 
Y 


l'algorithme, soit ici : 


NEM 
ios s (E Eu latino 

à ux)- a Bo ux) 0 
I)e (x)= a L(x)- ale d Q-l q= he Ge de 
i()- 15 ae p) ele; 1 y) p()- 261-7) 


2 
X X 


S(x) R'(x) e R(x)- e, (e 1e +1-7)-x (e - ae - B) 


2 
X X 


py) RG)u (x) S()u o) - -a p s()- 3a p (x (e -a e, - f)- (e, Ife «1-y) 
PDT Heul) ee De a fe, p) e(e 1 y) EL 1) (x (e, ae, p) ele +l y) 


X 


LAS (I ose) - 5 lSCenu e) Re) 


= rt -aJe - pisa e p) clin - BXa(y -1)-«) ylab aalders 


p(x)-- : 
etl 


1 


Ainsi les auteurs J.Letessier, G.Valent, J.Wimp prouvent que la fonction eer [eta -y4l 
S ZSNS 


respecte l'équation différentielle suivante (formule 3.21, x est changé en z et le tout est multiplié 
2 
par x//ei) : 


- i-e Be): - 4) 28), 


dz 
(a e ff e\a+p 1)22 «eo py e) (j Fla b+ e?)+a g)z«( +) .- ze. 


sa (a -eY8-e)z te NF -e)pe)- 


"nad rA 


a paz, AtB -e)- C to B jen B) Ge (e -7) 1o) 
( di (a-e8-«) m 


+ 


(e +I? -«) | 
a p| z+ LL |y(z)=0 
| (a-e XB -«) 
Par identification des termes en dérivée seconde, premiére et fonction avec l'équation différentielle 


zz E ee 1z a)+6 z(z a)+ez(z Neeg: q)y(z)-0 


2 T 
Z 


de Heun, 


2 


z(z-1\z y 8 +6+e)z°-(y(l+a)+6 a+e)z+y eE B z-4)y(2)-0 
a e-a) = ale +1-7) 
dz? Q e 8 «) (a af a) 


a p (e 1e -7) a B (e +e +1-7) 
(a-e p-e) (a-e p-e) 


Terme 


Terme ylz)> q= 


at) spen s EE Peg ce) ehe pe) B): 0e (s -7) 
Tene wu (a ni Ir -=ê 
Y Z 
on obtient : =(y+8+e)z -(y(1+t)+8+e)z+ya 
y*óte-zatf«l 


_F+)a(a-7)_rala+t-r) | ale +1-7) 

(a a\B e) (a a KB «) (a -eXB - e) 
e(a B1 -e)-(y «1 a Be’ ap. 
(a-e, (p-e) 


La dernière expression est automatiquement vérifiée à l'aide des expressions précédentes comme 
suit : 


ya 


yY(1*a)-8 a-e 


e(a* «Q7 -e)-(9 «1a Be? ea p yea) a. 


£--]l 
(a - 8-6) 
" ee +1-y{a+B+2) (a 2 2 y) es +1 y\a+B+2) 
eue) ONCE (x-e9-«) 


2 ea «107 -«)-(5 +1 a B*ej ea B=(a e\B-efi-y)-e(e *1-y Ya B 2) 
Terme ef >- *1)-(a* B41)-(1-y)-(a* B2) P 412 y 

Terme e/2y-y-l 
Terme e| 22y(a* B*1)- -(.-y Ka B)-(L-y (a - B 2) - 27 (a B 1) 


Il vient donc finalement un résultat assez simple : 


e. 8 £ La, MR 


dz? 2 EUR dz z(z -1z -a 
&--]1 
Contrainte de Fuchs y+0+e=a+p+l 
ee, -1— y) (e, * 1e -1— y) e +1 
q-ap =a pa 
(e, ae, -B) (e, -aXe - B) € 


dont la substitution des paramètres données par J.Letessier, G.Valent J.Wimp est la suivante : 


2 a Bz+— 
qus. ON C NEE: dy) y(z)=0 
dz z 2-1 , 1 dz 1 1 
z pg z(z- Ze 
1 S e +1 | 
a=t E q E s p F Contrainte de Fuchs y+8+e=a+p+1 
1 


C'est exactement la formule de factorisation de l'opérateur différentiel de la fonction 3F2 (avec les 
paramétres de J.Letessier, G.Valent J.Wimp) donnée par  K.Takemura (formules 1.9 et 1.10, 
attention une erreur de signe sur le paramètre q), sous la forme : 


e d? d = 
Lagya =55+ d + + zd + Bu < Opérateur différentiel de Heun 
UU gz z z-l z-ajdz z(z -1Yz - a) 


&£--] y+0+e=a+p+1sy+ô=a+p+2 
e (e, -1- y) q-a pte *1-Y) V pat se 12 $e e 8)+7 us qd €& 


(e, - a Ke, - 8) (e, - a Xe, - B) e l-a ape+l 
Factorisation  L,5,,,,, = E " e +1 i 1 - y 1 Hs, . 
z Z z-l z-a d 


C'est aussi la méme formulation donnée par T.A.Ishkhanyan, T.A.Shahverdyan, A.M.Ishkhanyan en 
formules 33, 34 et 35 dans laquelle on peut également interpréter la valeur caractéristique q de 
l'équation de Heun comme étant racines d'une équation polynomiale du second degré (formule 


32): 


S e +1 n ele -1— y) S aa p 
(E Ga ra aa 
ô=a+pB+2-7y 
LUE, Ceres) ve Butesiibey ls adea bea E, 
q-aaf 
2 ee *1- y) 
(e, - a Xe, - B) 


Si l'on observe bien, il s'agit de la méme équation du second degré obtenue pour les solutions de 
type (II.3), avec N-1 et en substituant t<->a : 

ô=a+B+2-y 

q racinede q^ *q(a-1-a(8-24y(20 -3))*ay(6 -2-a -a(8 -1f1« 7))- 0 


Autrement dit les solutions de type (II.3) exposées par T.A.lshkhanyan, T.A.Shahverdyan, 
A.M.Ishkhanyan ne serait ni plus ni moins que des solutions particuliéres de l'équation de Heun 
déjà observées auparavant, sous la forme de fonctions généralisées hypergéométriques, tout du 
moins pour les cas N=0 et N=1. 


L'identification des solutions de type (II.3) avec les fonctions généralisées hypergéométriques 
correspondantes est facile avec N=0 et N-1. Puisque pour chacune de ces solutions il s'agit des 
développements de Fróbenius d'indice O autour de la singularité z=0. Les deux solutions coincident 
donc à une constante multiplicative prés. 


Pour N=0, la correspondance est immédiate, pour N=1, on écrira la correspondance avec diverses 
relations entre les paramétres : 


Lgs mosso deals) CE aur 
ylzF,F la, B: L E (a) -1) F(a, fy; ) 


E ,acaa B+all-6) 
MN 
q racine de q-q(I-y*a(a* «2a B)}+aa B(I-a 1» B)- y)-0 
e (q-aa fi«ai- 5) a-aa B«(a-1 -y))- i-a la -rX1 p-7)-0 


re +l:7,e D 


|. acp B etl vs ee *1- y) d- (cal B)ey 
aa MEME T A or eer 38 


Qu'en est-il maintenant pour N-2 ? 


Je peux utiliser une méthode similaire à celle employée dans l'article « Some differential equations 
satisfied by hypergeometric functions » de J.Letessier, G.Valent, J.Wimp, soit calculer l'équation 
différentielle du second degré de la fonction hypergéométrique Æ (a, Be * Le *1;y,6,6, ; zy à l'aide 


de l'algorithme exposé par les auteurs, puis d'en déduire les relations entre les divers paramétres 
par identification avec les termes de l'équation de Heun. Toutefois pour l'instant je laisse de cóté 
cette méthode car le calcul de l'équation différentielle du second degré est suffisamment complexe 
(méme réalisé avec Mathematica). 


Je préfére la méthode de T.A.Ishkhanyan, T.A.Shahverdyan, A.M.Ishkhanyan, consistant à partir des 
termes du développement de la fonction généralisée hypergéométrique E, (a, Be +1l,e,+1;y,e,6;; z) 


qui suivent une récurrence à deux termes et de les identifier à la récurrence à trois termes du 
développement classique de Fröbenius autour de z=0. 


N A; =(j-1+a\j-1+8) 4, =(j-1+a\j-1+ p) 
X)-Xez ÀB=-jGlari)+e(-a)+ala+ph+y-D) | e7-N218, iG) (a7) Nan B)r-)-a 
j-0 C, «a(j Ay * j) C,=a(j+1fr+ j) 


Aic; * B, c, t Cie, EU 


Nous savons par ailleurs que par identification de la solution, le développement de Frobenius 
autour de z=0 est également la solution du développement (II.3) à une constante multiplicative 


prés. 
Comme les termes du développement de la fonction hypergéométrique „F(a, B. +1e,+1:7,e,e, ;z) 
suivant la récurrence à deux termes suivantes : 

rlari) ro) (a*je*)1 ., 6 1, (j-M8*j-U.— (a+ile+i) 


C.= >h, x 
FO Tap) Tti) ee 3 ^ e, j — (vei) — (e+j-Ke+j-1) 


En injectant dans la récurrence de Frobenius autour de z=0 d'indice O, il vient : 


Aj eo j AB j-1) 
(s À; 
N=2=4B,=-j(j(1+a)+2(a-1)+a(a+B)+y-1) q=+8,+Ch.=0 
C,=a(j+1\y+ j) | 
(*i-) eti-le+i-),, o 1 eriX8+i) fe jj) 
n . x ] jo Sept ] x ] 
atj-IÂB+j-1) — (aja +j) jio (wj) (e * je * j) 
e jx(y* j-1)x(e + j -1Xe, +j 1)+ B (e + j fe, + j)* ax(a  j KB * j)x(e, + j+1Xe, + j+1)=0 
e ji j-DMe + je, + j-1)+ Be + je, + j)+a (a+ JB j)xle + jte +j+1)=0 
Terme en j* = 1-(1+4)-a=0 
Terme en j 2 39e +e, *y *3-e-e, a(2-*e, * e, * a. 4 f) y*a(2*e *e, +a+B)=0 


SUM 


C'est une équation sur l'indice j à l'ordre 2, en effet le terme en j* et le terme en f? s'annule 


directement. 


Il reste donc un système de 3 équations algébriques sur les paramètres a, q, a, 6, y, el et e2. La 
première équation algébrique correspondant à l'annulation du terme de degré 0 en puissance de j 
donne la valeur : 

(e, «1e, +1) 


ee 


q=aaß 


La deuxième équation algébrique correspondant à l'annulation du terme de degré 1 en puissance 
de j, donne la valeur : 

2 ee, *(y-1ft-« -e,) 
a p (e(1+e )+ eee ))+ ee, ee, -(a * BL e, +e,)) 


d - €; 


La troisiéme équation algébrique correspondant à l'annulation du terme de degré 2 en puissance 
de j, donne la valeur : 

3+e +e, -2y 
a p(l+e +e, )+ee le +e, -1-2(a +B)) 


d - €; 


Toutes ces expressions sont à comparer avec celles données en formules (40) et (41) dans l'article 
« Expansions of the Solutions of the General Heun Equation Governed by Two-Term Recurrence » 
des auteurs T.A.lshkhanyan, T.A.Shahverdyan, A.M.Ishkhanyan. L'expression de a est différente, 
mais cela n'indique pas pour autant qu'elle soit fausse. 


Il semble donc y avoir surabondance de solutions pour le paramètre a. Qu'en est-il en réalité ? 


Toutes ces expressions sont intéressantes en ceci qu'elles permettent de définir les paramétre a et q 
en fonction de la somme e.+e, et du produit eje; comme ceci : 


q-aa f l+ py tsp 
5n 
345, -2y 
a Bess )* ps -1- Xa + B) 
2 p; *(y- Ml -s;) 
a P (sall +s9)-2P2)*+ Pol Pa -(a * BXLe s) 


$,76 t6 pyc-e£—38-py, 


a= py 


Les deux premiéres équations s'inversent facilement pour donner : 


_3+a-g+a(a B+2(a+$6) 2+a(a B+2(1+a+8))-27-q 
"E a-l PEAU (a-1Kq-a a f) 


L'injection de ces deux expressions dans la troisiéme expression permet de factoriser le polynóme 
en q suivant : 


P(a)- (aaa BJ +(q-a a B)(4a-a(3+a+B)-2+y)+2a(a-1a Blq-a a B-2a(1+a+p)-2+2y}+2a(a-1fq-a a BXL- a) B) 


Ce qui est exactement le polynôme dont q doit être la racine pour assurer une développement de 
type (111.2) fini. Cela signifie que la troisième équation algébrique est automatiquement vérifiée par 
la factorisation de ce polynôme. L'expression (41) de a : 


2ee, +(2-yÎ1+e +e,) 
a plite} +(1+e,) -1}+ee,(2 ee,-4-(3+e +e fa ^ B -1)) 


a=ee, 


donnée par T.A.lshkhanyan, T.A.Shahverdyan, A.M.Ishkhanyan est différente mais en l'injectant 
également comme autre équation on factorise encore le polynóme Pq): Cette expression est donc 


juste. 


Poussons toujours plus loin (vers l'infini et au delà ...) avec N=3. Pour la détermination de la 
fonction généralisée hypergéométrique suivante: y; F (a,p,e «Le cle +1;y,e,e,e;z) dans la 


récurrence de Fröbenius d'indice O autour de z=0, il vient l'identification suivante : 


jo S l eri-N8+i-t), la+ jeje +) 
a À ti) (eeiclke rites j-1) 
A, - (a j- 1 j-1) 
fs A; 
N23 et 1B, 2-j(j (Lea) N (a-1)a(a B) y -1) q> +B EC -0 
C, za(j Xr * j) 


j 


On obtient une équation algébrique de degré N=3 en j, qui par annulation de chaque terme de 
puissance de j donne un système de N+1=4 équations algébriques, dont celle de puissance O (la 
première équation) donne invariablement l'expression : 

(e, +1{e, +1e +1) 


CCS) 


q=aa p 


Toutes les autres équations algébriques donnent une détermination de a en fonction des 
paramètres e,,e,e sous la forme a-f(ee.e) j=1.,N=3 Ces N=3 équations permettent en théorie 


de déterminer les trois paramètres e:e,e; en fonction des paramètres a,q,a,6,y, N=3. Il faut plus 
précisément tenir compte de la remarque suivante : 


Remarque importante : étant donnée que la donnée des paramètres ee- et e; est indifférente à la 
permutation quelconque de ces trois paramètres, il s'ensuit que chaque expression de a est 
également fonction des paramètres symétriques s;,s; et ss suivant : 


S 56 te t6 
S, = €€, t 66, t €€, 


53 606 


On obtient alors un systéme d'équations algébriques sous la forme : 


]ts +s, +s 
q=aaß l 2 3 


3 

s,(6+5 -3y) 

a = 

ss, +a B(l+s,+s,)-3s,(1+a +8) 
s(1-s,+5,+3s,-y(1-5.+5,)) 
35 -s(a+B1+s,+s,)+a Bls,+55, +5, -35,-255,) 

s, 5,+25,-4-7 Q2 -3) 

s (2s, - s, -1- (a B (25, +3))+a p Is +55, +57 -35.) 


Les trois dernières équations permettent d'obtenir les expressions de 5:,52,53 suivantes : 


_6-qg+3a(l+a+p}+aaB-37 
a-l 
q a(a * f 2a f) y): 
«qT«a(p (4+38)+a (5+88+48°)}+2a(2+98+48°)+107-37°+2a(-1+a+B-2a B (-1+y)+7) 
1 -aa Baa +a°(3+78+26°}+a(13+188+78°)+19y-37" -a (C1 18a (158 C17) ey) 
"ir gel q «aa B(a( «a^ Ba p) 2 y) qli+a(a+B+2a f) y) 
_aaB 2«a(a B+2(1+a+8B)) 2y-q 
a-l q «aa Blal3+a+B+a p) 2 y) qli+a(a+B+2a B) y) 


1 


s3 


La quatrième équation est en surabondance pour la détermination de s,,555; en fonction de q et a, 
mais elle est automatiquement vérifiée grâce à la factorisation du polynôme P:(q) dont q est la 
racine. Le polynôme P:(q) du développement de type (II.3) est le suivant : 


qg-2q(5+a(2+3(a+B)+2a B)-3y)+ 
a(4+18B+11 f^ «o (18 B+6 B )+2 a (+22 8+9 pP) 
+q/1-2a(-6-158+4y7+118y+a(-15-10B+11y+98y))+ + 
+33-40y7+117° 


a 6838 P )ea «pop +28" ea (^ d 38? 4 vuU. Enea 


P(g)-1«2q| «a'(o (9 p (83-417)+ p(15-227)-9 y «a? (0 p (15-22) (5-9 7)-97)-98 (1+ B(-1«7)-3C 6115-67 y): 
a((a- Bo-5y y a B(2-25y 117) 

a [611a 62 ea 61186? p?) 
22? (18-9 (o? + «oTt af Ya B)-37 a B+30 p? ]y 
«ay (18(1+a+BXy-1)+a B (11 -10) 

-6yQ-3«r) 


b 


Par extension pour n'importe quelle valeur de N, il vient l'algorithme suivant de construction des 
paramètres de la fonction généralisée hypergéométrique :y(2} F. (o, ,e «1... 1i y ese, iz) 


solution de l'équation de Heun : 


INNEN! 


lon À — (yj) I +5-1) 


A; (a j AB j-1) 

TE À; 
B,- dj Ga) N (a 1)+a(a+B)+y 1) HE OMAN 
C, -a(i My j) | 


On obtient une équation algébrique de degré N en j, qui par annulation de chaque terme de 


puissance de j donne un système de N+1 équations algébriques, dont celle de puissance 0 donne 
I-N 


[It +1) 


invariablement l'expression : q=aaß TAER Toutes les autres équations algébriques 2 à N donnent 
IE 
l=1 
une détermination de a en fonction des paramètres S,Sz...,Sn sous la forme a= f(s.s,...s,) j-L., N- 
j , LS Pos, geras 


En effet toutes les équations algébriques obtenues sont parfaitement symétriques à toute 
permutation des paramètres eye; , .., en il s'ensuit que chaque équation est fonction des 
paramétres symétriques : 


I-N i j=N i j=N 
$, 71 sl s= $ee, s= X eee 
l=1 ij-li»j ij k-Li» j»k 


is IN 
Sy 5 Ys X..X6, Sy =| [e 


isiy => ly l=1 


L'introduction des paramètres symétriques S,S2, .., sy permet de linéariser le système initial 
d'équations algébriques. En effet les produits symétriques contenus dans les coefficients hj; 
deviennent des sommes : 


I=N IN 
etj Sy j! 
„ Lue i a i). [+2 EN TEE 
NET ; I-N , ; IN 
j ytj-l ; jj y*tj-l 

I-N I-N ; I-N I-N ; l ) [1 et] -I) ( ) x -1) 
| (+5) Ys Merret 6-9 > o | 
7 i i (a+ JB jx} syal +1) 

HT Fe j 


(j + Wy + nr Sy4J 


1-0 


D'où une équation algébrique de degré N42 enj : 


IN 


jtsi- Dsl- lat IB iZ sr GH) 
A b kei 1)«a(a* B)*y-1)-q*a y =0 
Naya Zs 
1-0 
js j- DA G (La) N (a- 1)*a(a* B)«y -1) +q) D j'+a(a+ jKB- j)x xs) £l) EU 


1-0 1-0 


On vérifie aisément que le degré N42 et N+1 de cette équation algébrique s'annulent : 


l=N l=N 


+ sl- -U U (l+a)+N (a-1)+a (a+ p)+y-1)+q)x 2 +a(a+ j ^ j) DUE j+1) 20 


1-0 


Terme j^ -s,-(I-a)s as, =0 


Terme jug ERI N)+5,(y 1) (1-2),  s(N (a 71) a (a B) y -1)) a (s, 5, s(N -1-a * B)) 
=s -(I*a)s +as 20 


L'équation algébrique en j est donc en réalité bien une équation de degré N en j. L'annulation de 
tous les termes de puissance de j aboutit immanquablement à un système de N+1 équations 
linéaires dans les variables 5:,5:, ..., Sn. Les N premières équations servent à déterminer 51,52, ..., Sn 
en fonction des paramètres a,q,a,6,y,N : sz gl frag N) j=1.,N- La N+1 équation linéaire 


détermine également une valeur ,- fils...) S dans cette derniére expression on injecte les 
valeurs s,=g;(a,8,7,a,9,N) j=l, N? alors on parvient à factoriser le polynôme caractéristique du 


développement de type (II.3). Cela permet de vérifier que la N+1ème équation algébrique est 
automatiquement vérifiée par la condition que q est racine du polynóme caractéristique. Ce dernier 
polynôme d'ordre N+1 en q est déterminé par l'annulation du déterminant suivant : 


[Neg SEU. A "us 5 oue 0 

a Laer N- DEAN Ny) Au. cedro dione cos 0 
j y-N+j-1 0 A3) B" -q cu» 0 m 0 
BÍ aa Be aif Ne j-l)ra(N- iere) De - oom 
cuz -(a- If, - N^ jXj«1) B $5 E E S P 0 
» " 0 40» puis ci 

| 0 a ioci. XR Cup Ru 


Remarque : en repartant de la récurrence originale sans poser ; __y 


4; =(j-1+a\j-1+8) 
B,=-j(j(a+1)+e(1-a)+a(a+B)+7-1) 
C,-a(i 1f j) 


Alors l'équation algébrique de degré N+2 en j s'écrit : 


itrsi- less j+a)-e(a-t)+ala+p}+y-1)+9 Ys talat IB xs H1) 20 


1-0 1-0 


dj 


Et c'est l'annulation du terme en N+1 qui permet de déduire la contrainte :, .. X 


l=N l=N 


IN 
i+j-È sG- -GGea)-e (a-1)+a (a+ p)+y-1)+4)x} 5, +a (a+ IB Ts, G1) -0 
1=0 Io 1=0 
Terme j"?=s,-(1+a)s,+as, =0 
Terme re s, sy N)+s,(y 1) (1-2), (76 (a-1)*a (a * p)+y-1))+a (s, s, s(N -1-a ^ ]) 
=s -(I*a)s, +as +s,(N+e)(a-1)=s (N+e)(a-1)=0=e=-N 


Il reste une toute petite étape qui consiste à revenir des produits symétriques aux variables e; . Il 
suffit pour cela d'utiliser la formule du polynôme de Vieta : 


j=N j=N 


%=1,5; j=L.N 2 Vietalt, N)= Ÿ 53. ( 1)" #1 = ( e) 


150 j=l 


Par un programme numérique, il suffit donc de rechercher les N racines du polynôme de Vieta, 
calculé à partir des valeurs s; pour obtenir toutes les N valeurs des variables e; . A titre d'exemple 


voici des valeurs particuliéres des paramétres : 


a £--N 
er ET 


s 10411 , 453403221 , 125887503663 , 1625509485739322]  176651269997811147 
OS uS" s4000 7 zu9D80 C * 1160858560000 IBTI7736960000 
)=0 & q- (23.3539,14.4862 , 7.81156 , 3.32818 , 1.07301} 

s,=101.527 fe, 258.8103 

s, =2898.81 |e, 29.9684 

s, =228729 © Je, =123811 

s, =8019.87 le, =0.367528 


6 7 3 7 
5 


Pour q=3.32818 > 


Et les graphes respectifs des fonctions construites coïncident parfaitement : 
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0.99 HeunishkhanianExpansion(nTerm,a,g,a,8,y,6,€,z) 


HeunlshkhanianExpansion(nTerm,a,q,a,8, y,6,e,0) 


—— gFs(a, B, 1«e1, 1+e2, 1+e3, 1 e4; y, e1, e2, e3, e4; z) 


Il se peut que les paramètres e; racines du polynôme de Vieta, soient parfois complexes. Dans ce 
cas, le conjugué complexe e; est également une racine du polynóme de Vieta. Dans ces conditions 
la fonction généralisée hypergéométrique reste à valeur réelle. Partant de la définition de la 
fonction généralisée hypergéométrique, et considérant par exemple e1 et e2 comme les racines 
complexes conjuguées, il vient : 


EN PUR DRE ACTEUR T() Wn Trey). Tle) 2 
| i Y, eize €x T(i+e,) Tle, +1) I! 


1-0 T(a)r(8) T(y«1) ja 
rer) ind E 


a, BA el te. e, er(a«n)n(gs)) r(y) Fee eren +1) 
Es à z[- XL. si a x 
AAC & Tar)  r(y«i) (te rll«a) T(e, «1)rle +1) j3 rte.) Te, +1 


24N X 


= 1*e,l1*e,,..] e, 
Comme ree el) (ne +1 x^ prp ral à valeur réelle 
/, 6€ Ex 


L'existence d'un tel résultat n'est donc plus une conjecture puisque le systéme d'équations a été 
linéarisé. L'algorithme est applicable sous Mathematica et donne toute satisfaction, seulement 
limité par le temps de calcul (qui s'avére négligeable au final, si le probléme est totalement 
linéarisé). 


Déduction des solutions pour & entier positif 


Un article de 2020 du méme auteur A.M.Ishkhanyan « Generalized hypergeometric solutions of the 
Heun equation » apporte un complément de deux natures. Tout d'abord en plus de l'article de 2018 
des auteurs «Gauss-hypergeometric expansions of the general Heun functions governed by two- 
term recurrence relations » dont on a déjà parlé, | A.M.Ishkhanyan introduit les produits 
symétriques s; et permet donc la linéarisation du systéme d'équations algébriques et donc la 
solution compléte telle qu'on a pu le voir. 


Ensuite il traite le cas des & entiers positifs en page 6 et 7. Pour cela rappelons une des 
transformation homotopique de l'équation de Heun qui permet de « renverser » le signe du 
paramétre & : 


y(z)* (z- a) * yu (z)-9 »y(z)= HeunG,(a,q y( £);a +1-e,B+1 £,y,0,2—8;2) 


Alors le schéma de construction est le suivant : 


Jeux de paramétres Racine qu et ee, 
E 
atl-e,B+1-6,7,0,2-8=-N ' |yg(zey Fol 1-6, f 1-6,0) 1.6, 1; Y, 6,6, ; z) 


e y(z)- (z-a) ™ , ,F, (a -1- N, B-1- N,e *1..,ey 1; y, 6,6 ;z) 
Solution de l'équation de Heun avec q= qq +y(1+N) et a,ßp,y,ô et e=N+2 N20 


Propriété supplémentaire du développement de type (II.3) : la singularité « apparente », disparition 


du terme logarithmique du développement de Fróbenius autour du point singulier régulier z=a 


Un classique résultat de la théorie des équations différentielles ordinaires de degré quelconque 
stipule que lorsque l'équation indicielle du développement de Fróbenius autour d'une singularité 
réguliére zo donne deux racines dont l'écart est un entier alors la seconde solution indépendante 
autour de la singularité zo; comporte un terme logarithmique. Il existe des conditions biens 
spécifiques pour que ce terme logarithmique disparaisse. Lorsque ces conditions sont remplies on 
dit alors que la singularité autour du point régulier zo est « apparente ». 


Pour le cas de l'équation de Heun je m'en référe pour partie au formalisme de la théorie de 
Fróbenius (exposé dans et pour partie à des considérations spécifiques portant sur les 
développements de Fróbenius. Plus précisément considérons l'équation de Heun et les diverses 
exposants de l'équation indicielle : 


ves(L-2 ye i-o 


-] z-a z(z-1Yz-a 


Condition de Fuchs y+8+e=a+fB+1 


Point singulier z=0 p=% qy-0 I(r) r ?er(y- )or 


Il 
© 
x 

Il 
eR 

l 

ES 


Point singulier z -1 p= qo 02 1(r)= r ? e r(5 — 1)=r=0 r-l1-ó 


Point singulier z=a py-6 Qo 02 1(r)= r ? e r(g— 1)=r=0 r-l-& 


Point singulier z =% = r^ +r(1 -p)*4 = (r- ar- B)- 0rza rz 


Et regardons spécifiquement la condition pour laquelle le terme logarithmique sera absent, pour 
un développement autour du point singulier z-a,. Étant donné qu'autour de la singularité z-a 
l'écart en valeur absolue entre les deux racines indicielle est de | 1-£|, il suffit que € soit un entier 
négatif égal à -N pour que le terme logarithmique soit potentiellement présent. Les deux racines 
soit alors respectivement : rı=1+N et r,=0. Les fonctions f; (voir théorie de Fróbenius exposée par 
E.L.Ince dans son ouvrage « Ordinary Differential Equation », pages 356 et suivantes, page 396 et 
suivantes « Solution of linear differential equations in series » et paragraphe 16.33 page 404 et 
suivantes) sont définies comme suit : 


Ó- [L2 o) x Cin 


z mm z(z -1) 
(z a)P( JE H à 1) (y(a 1)” H ZU 


f(z-a,r)» r(r 14 g) z—a ( yer aXa ne (6 23 apja ed; 


— fr)» r(r -18) 


[| ja (r Ya 1) «(ró -af)a' " 
se Cp ere Dr eser. uus 


Lorsque £--N, r;-1*N et r,=0, la condition d'annulation donnée par E.L.Ince du terme logarithmique 
s'écrit sous la forme : 


AN) f(N-1 f4(N-2) . f) Soal) 
AN) f(N-1 f(N-2) . f) flO) 
0  A(N-1 f(N-2) . . 


Fa =n =0)= 0 £0 (0) » 
AQ) AO) 
| 0 0 o A) (0| 
fN)- -N 
Avec Fe TE Worte ee ar eple s 


Lorsque e=0 alors F(0)= /(0)2 0 & q(a-1)+(-q+aBa=0& q=a ap . 


[4 50). i i 
Lorsque £=1 alors ? LA Po AOAO AAA (0)- 0 


expression donne une équation du second degré en q : 


. Le numérateur simplifié de cette 


q'-q(-y*a(y*ó-2«2aB))* aof(a(y - 8 -1-oB)- y) 20 
Étant donné que la relation de Fuchs a court : Y +ô-2=a+ B, Cette équation devient : 
q!-q(- y -a(a B*2af))* a aB(a (1- a X1 8)- y)-0 
AQ) f) AO) 
fi 


Pour e=2, alors (r-0)- f) fii) f(0)|-o 
o 40 (0) 


il s'agit de trouver les racines d'un polynóme du troisiéme degré en q, sachant que la relation de 
Fuchs a court : Y *9 3-0 * B. L'expression de polynôme en q’ est la suivante : 


a^ GB(- B)-2a(1-58 «38? )« o^ (0-68 3p?) 
q'-q'y -4-a(1-3(a - B a) q| +a(a(4+48 -4y -6p y )- AB(y -1))+ 

+2(y -1y -2) 
-a apla’ (a? -3a - 218? «38 -2)- ay (4 - 4a 4B +3aB)+2y(y -1))=0 


Cela coincide avec l'expression (2.7) obtenue par K.Takemura dans l'article de 2022 « Heun's 
equation, generalized hypergeometric function and exceptional Jacobi polynomial ». Cette 
expression peut aussi s'écrire sous la forme (voir article de 2018 T.A.Ishkhanyan, T.A.Shahverdyan, 
A.M.Ishkhanyan « Expansions of the Solutions of the General Heun Equation Governed by Two- 
Term Recurrence ») : 


(la-a apy +(q-a afi (4a - 2- a(3*- a + ß)+ y)*2a(a -1of fa -a af —2a(1+ a B) 2(y -1)) 
*2a(a -1(q -a ap X1- a X1 B) 0 


Les deux développements indépendants pour les deux racines indicielles r1=1+N et r2-0 de 
l'équation de Heun autour de la singularité z-a (ici apparente) sont théroiquement construits 
comme suit : 


Zeiex)-G-ay" Seed) — e De cr) EOD 
3 [Ie 


Zen) 


=D eO), 6-2) 


La relation de récurrence pour les coefficients du développement est la suivante : 
jeN g(r) rv) ga (r)firv —1)+...+ a (r)f,(r+1)+ Cf;(r)= 0 


Cette expression de la relation de récurrence conduit à un développement de la solution 
suffisamment compliqué pour ne pas poursuivre le calcul plus avant. Mais on peut simplifier la 
récurrence en exprimant l'équation différentielle sous la forme suivante : 


acc e sa ect af? m presen: pcs 


-1 
Soit p, (z)y"(z)+ Pi (z)y'(z)+ P> (z)v(z) =0 où p, (z), Pi (z) P3 (z) analytique autour de z = a 


z jato . M . P4 M . 
Le développement de la forme : Ga) ce, (z-a conduit à une relation de récurrence à trois 
J 
j=0 


termes, comme suit : 


P,=(r+a-1+j\r+B-1+ j) 


up Q, --(r jy -q-(r+ jXr-Me y +e)+al2(r+ 5 ta (r4 ja +B-1+6) 
R, =aļla-1\r+1+jÁr+e+j) 


Pc;,+0Q,;c;+R;c 


Au passage l'équation indicielle se retrouve en annulant le degré minimal soit celui du terme (z-a)' : 
Terme (z-ay -0e P, -a(a-1)r(r-e41) 20e r(r+e-1)=0 


D'autre part pour la racine indicielle r=0, le coefficient R; dans la récurrence Pe. (5 Qc - Rc, -0 
HIT PC: vi 


s'annule lorsque j--&. Pour que le coefficient c;.; demeure fini et qu'ainsi le développement d'une 


j+ 


solution de la forme S c(z-a) puisse exister, il convient que l'application de la relation de 
J 
j=0 
récurrence pour j=-E=N (lorsque £ est un entier l'écart des racines indicielles est alors un entier) soit 
vérifiée : 
€ -—N P,c,,+0O,c,+R -0 
P,=(a-1-eXB-1+e) Q,=aaB-g+e(s+1{2a-1)+e(y+e-a(a+B+1+e)  R,=0 


(a-1-EsXB-1+e)c.,+(aap -q+e(s 1 (2a 1)4 ely HE ala + B+1+ ec =0 


-eĈ-e+1 


Le terme <- est donc parfaitement défini comme suit : 


T" -(a-1-eXB —1-2) B 
* (aaB-qg+e(s+12a-1)+e(y+e-a(a+B+1+e)) ^ 


ainsi que tous les autres termes du développement , j >o en appliquant pour ces derniers la 
Ln 

relation de récurrence : 

DuC kA + (OTE R z 0 


—€e+j+1 F R —E+j 
—£-ctj 


P,,,c +O oaj ,,;+R,,,c =0c 


—e£+j“-e+j-l —eE+j“-e+j+1 


Si par construction les deux types de développement eS d(z-a) et 
j-0 


Î=+0 


y, (z)= Y c (z-a sont indépendants alors la condition précédente garantit qu'il y a absence 
J 
j-0 
complète du terme logarithmique dans la solution indicielle r=0 (racine la plus petite) autour de 
z-g. 


Or une solution autour de z-0 peut étre est construite sous la forme d'un développement fini de 
fonctions hypergéométriques 2F1. Lorsque z est réel et si le paramétre a « 1 alors la fonction 2F1 
est bien défini par son développement classique. Si la paramétre a » 1, il convient d'utiliser la 
continuation analytique de la fonction hypergéométrique 2F1 (utilisation de sa représentation 
intégrale par exemple) pour s'assurer de l'existence du développement de type (II.3) également 
autour de z-a. 


Comme cette solution est indépendante de la solution indicielle r=1-e, alors là encore si l'on 
retombe sur la méme condition polynomiale conduisant à des valeurs spécifiques de q racines de ce 
polynôme, on peut affirmer que l'annulation du polynôme caractéristique Py(q) est bien la 
condition d'absence du terme logarithmique. 


Je me suis assurée que les polynómes Py,(q) calculés soit par la condition de E.L.Ince, soit par celle 
de K.Takemura et par celle donnée par T.A.Ishkhanyan, T.A.Shahverdyan, A.M.Ishkhanyan sont 
rigoureusement identiques tant que les valeurs de N sont compatibles avec la durée du temps de 
calcul réalisé par Mathematica. 


Un cas particulier de récurrence à deux termes lié au développement (11.3) 
. H j=% 
Parmi le développement Type (11) (= Ye, .Fifa Bira c) 
j=0 


dont la relation de récurrence à trois termes sur les coefficients c; est la suivante : 


-0 


(gt) q)c, ' cHe 
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Type (I) A, 


40 = a(a-y,- j:M(B-y,-j« (y - v, - i1) 
: yy tj-l 


B" =aa f (a Wy, 7j y «c Yo j)*a(y Yo jfa*p Vo j) 
C z(a-1 y, e jy £77 i71) 


ju 


d : 7+0=2 q-aa «al 8)e Vo=7+e a+B+1=y+0+e 


si B” -q=aa B-[a-Iyy, jy 67-529 -n- Ya 877.) 
TOUS posons = pun. o af je j(a- iy 127 j)-a(e jX9-1- j) 
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La récurrence de type (II) devient une récurrence à deux termes, et l'expression du développement 
de type Il devient : 


reca i) e- Bi) j) 
(y *e* jy e jelXj2) 


j: k=% 
yz)» Me, Fo, Biy +e+j;z)= Yes Ala, By +e +2k;z) 


j=0 k=0 
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ge. "un prn „Fla, sy * e 2k;z) 


2 2 
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240) Qr) 
Cy, = A, Cy avec 4j 


La valeur en z=0 est la valeur de la fonction hypergéométrique généralisée en z=1 : 


ux ux i) p (=£) B 
y0)- Y 2 Pi ANEN Febry S IM 2 AN? k- F 
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Développement d'une récurrence à deux termes à partir du développement (II.3) 


Dans l'article « Expansions of the Solutions of the General Heun Equation Governed by Two-Term 
Recurrence » les auteurs T.A.Ishkhanyan, T.A.Shahverdyan, A.M.Ishkhanyan suggérent un autre 
développement récurrent à deux termes afin de le comparer au développement 2.3 (voir formules 
(51) et (52) de l'article). Pour cela il calque la récurrence à deux termes sur celle des coefficients du 
développement autour de z-0 de la fonction généralisée hypergéométrique suivante : 


vaoFyalf *£-0,y 4 &- B,e +l. ey oi y +E, epey) 


dont le rapport entre deux coefficients successifs est le suivant : 


eed en et h=” zi yte 2] re Btj 1), E 


0 Ca J (ye j-1) 


Et que l'on injecte dans la récurrence de type (II.3) sous la forme suivante : 


ans 26e y-€ j+1)(B y-€ jy d j) 

/ ytej-l 

B? =aa -q« j(a-My*ej-1)-a(e* jfa* B-y-e-j) 
cun (a Wy «e jj 1) 


pu» n CU =i 


j=% 
yz)» 3e; „F(a, B:7+e+ j;z) (13) 
j-0 E À; 
AU + BU) c, £C, =0 h 


J 


, pino ; CU, =0 pour j>0 


! 
Comme auparavant utilisons les produits symétriques s; des paramétres e; et leur linéarisation : 


I-N i,jN i,jN 


ace $7 Dee, $,— $ eee 


l=1 ijj i, j,k=l,i> j»k 


I-N IN I=N I=N 
> (o * j)2 Ys, j I(e*i-02Y5, ES 
isiy m l=1 1-0 l=1 1-0 
Sy = Ys X...X €, , Sy = € 


isiy => div l=1 


Les rapports de deux coefficients successifs donnent : 


l=N l=N 
(y&e-a j-Myes- B j-1)xY sj (y+e-a+ jy *s- B ix Ys, G1) 
h,= à io h = [n] 


J j+l I=N 


IN N 
jore jx sy (j -1) (il ye j)x Ys, uj 
1-0 


1-0 


D'où une équation algébrique de degré N42 en j pour tout indice j >=0 (méme l'équation pour j=0 
peut-être retrouvé en posant j=0) : 


1-0 
A09 
B"? ag p-q+jla-Wy+e+j-1)-a(e+jfa+p-y-e-j) + + CP =0 
IzN d 
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IN N 
2 -ja(e+j- Es 1] +B E s, ij M (1-afy ea jfy+e-B+ Ts G1) =0 
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E 1-0 l= 
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IN 


IN 
e-ja(e* j-MxY ss G1) «(aa B-q+jla-Wy+e+j-1)-ale+jla+8-7-e-xS sui +ll-afy+e-a+ ree B i) s 1) =0 
1=0 


1-0 


On vérifie aisément que le terme de degré N+2 de cette équation algébrique pour un indice j>0 
s'annule : 


IN 


-ja(esj- srli- taa B-q+ a Wr j-)-a (et ie 87r t D 


1-0 


IzN 
*(I-ayy *e-a jy 7 B j)x Y sy (1) =0 
1-0 
jaler j-i 1) s -0)") 
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( aXy*e a+ n Fe - fà xs sj) 
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(l aly HE a+ jYy te- fà ju Qs? +N j)esGeN-1) 


Terme j“? 2-as,*Qa-1)s, *(1-a)s, =0 


Terme j" -s(y*e* N41-a-f) 


Le terme de degré N+1 doit également s'annuler, il vient donc une contrainte directe sur le 
paramètre 6 : 


y+e+N+1-a-B=a+f-ô+N+2-a-B=-0+N+2=0=65=N+2 


Il reste donc un système de N+1 équations linéaires selon les variables symétriques s;, ..., Sn. L'une 
quelconque de ces équations est en trop (prenons celle du terme de degré N en j). Elle sera 
automatiquement vérifiée par une condition d'annulation d'un polynôme en q de degré N+1, une 
fois qu'on lui injecte les solutions s; , ..., Sn en fonction des paramètres a,q,a,6,y,N obtenue par 
résolution du système linéaire de N premières équation (degré 0 à N-1 de l'équation algébrique en 


j). 


Pour obtenir la condition sur le paramétre q, il suffit méme d'écrire le systéme d'équations linéaire 
selon les variables symétriques dans l'ordre inverse sy, Sn-1, ...,51So71: 
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Terme j'— nee -0 
1-0 
à IN 
Terme j 2) d,s,,=0 
Jo Dds Sa )-o 


IN 
Terme j“ > Y dy j$y., = 0 


1-0 


Seuls les éléments diagonaux de la matrice | ü | en question se présente sous la forme : 
Jl 


d, 7 f;(a,a,b,7,N)-q' Les autres éléments non diagonaux sont uniquement fonctions des paramètres 
j , , 34:79. 


a,0,6,y,N. Cela assure bien que l'annulation du déterminant donne une équation polynomiale de 
degré N+1 en q. Une fois déterminé les valeurs de sy, Sy, …,S1, la forme du développement ainsi 
construit est le suivant : 


C € (y+e-a+j-1{y+e-B+rj-1) 4" 
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l 
On peut montrer que quelque soit N, les valeur du produit : cu et de s; sont invariablement 


égales à : 


s-Yecceeelren) ND jp) 


I=N IzN 


(te) 1+5, 


i 0487. -q+(N+1)ay+(N+1-a\N+1-8) 
TI Sy (a-1\N+1-a\N+1-£) 


Remarque : je suis dubitatif sur la convergence rapide de cette forme de solution, car une fois 
testée numériquement la solution ne semble pas de premier abord respecter rétrospectivement 
l'équation différentielle de Heun. 


Mais un test numérique a finalement révélé que l'équation différentielle est bien respecté comme 
suit. Introduisons l'opérateur diffézentielle de Heun : 


Diem. s c) eee) (+ [e e reo LEEA 9 


z z-l z-a (z-1Xz-a) 


L'application de cet opérateur sur les fonctions de base du développement se présente sous la 
forme suivante : 


PHeun(j,-1) ,F(a, B;y - & * j-1;z)+ 
Dhel £, Fs Biennial + PHeun( j,0) „Fla, p;y+e+ j;z)+ 

+ PHeun(j,+1) ,F(@, By +e + j+1;z) 
PHeur( j,-1) - (1-a) j (y +e+ j-1)= CHP 
PHeun(j,0)= j’ (2a 1)+ j (1 H y (a 1) (2a l)e aô) q+aapß aell 8)- 80") 


SN. | ; OD | (13 
PH. +1 H ó -1 A: 
eun( j,- ) a (c | J ES jH 


Expression obtenue en comparant avec les termes de la récurrence : 


A9 =a (a TE j+1)(B y-€ ju E j) 
t yte j-l 


B'3)2agf q* j (a Wy ej 1) a(e+j\a+B-y-e-j) 


ce --(a-Ay «e jj 1) 


T» d 
Soit: Diten| ; Fe Bree) pe to Ala fiy o6 j-hz) e BU? Fa, Bytet jc) AU! Kla, fiy ee jhz) 


PE j- j j+ 
z 


Alors si l'opérateur différentielle est appliqué sur une approximation du développement avec M+1 
termes, alors il vient l'approximation suivante de l'équation différentielle de Heun : 


jM 


; l T(y*e-a j)T(y -e- B * j) T(y +e) Iz 
M = QE . . =] Zt 
Vnde po esr) alo SU Taea ree) P 


e 
l=1 
DH £ u i cy AP x Fla, Bye £M -1;z) ec, (BU?) x „Ela, B; y+e+M ;z)e AN) x F(a,B:y+e+M+1:2) 
: Hu z (z-1Xz-a) 


Autrement dit l'équation différentielle est respectée modulo un terme résiduel qui demeure non nul 
méme à très grande valeur de M et qui ne dépend fondamentalement que des deux derniers 
termes de développement à M+1 valeurs. Soit que la convergence est extrêmement lente. 


Les premiers développement pour N=0 et N-1 


Prenons N=0 


Le traitement du cas N=0 est spécifique car le système d'équations linéaire n'en comporte plus 
qu'une : 

Sachant que a * B. -y-£-1 e y4e-a 2 B-1oy*e-Bza-1 

pus cu. 2g h _(r+6-ay+e-8) 

0 0 i 1 y+e 

e laa B-q-as(a-^ B-y-c)) (Ql -aky *e-aYy - B) -0 

> q={1-afy+e-afy+e-B}+aa B-as(a* B-y-«c) 

eq-(f-1fa-1)-a(8-1(a-1)*aa B -ac & q2ay*(B-1fa-1) 
Donc le polynóme de degré 1 en q dont les racines déterminent les valeurs de q est le suivant: 


j= 


xe Xe silere) qzaysla-iyp-) e c trt I tepti 


a Ga j(y+e+j-1) 
] _1 T(ree-as j)(re- Bj). Tree) 1 (y+e-a) (y+e-B), 
> C ji Cj = X X : - X 
jh Ify+e-a)l{y+e-f) — (sj) j! (yc), 


Symbole de Pochhammer. (y *£-a), (y +e p), (7+ e); 


J 
Compte tenu de la relation de Fuchs on peut écrire encore plus simplement 


g=ay+(a-i\p-1) Relation de Fuchs a+B-1=7+e 
Symbole de Pochhammer | (a -1) ; (B 1), (a 1), 


X-F : Et e, „Fla, Bia+B-1+ j;z) 


mj! a+B-1), 
(a S LU 7D, 6-0, Ff eL Bea B ir) (our. x ae cce xd. Me 
5e) Aj (a+B-1), a+B-1+)j » re ee z-a | e rar per HA : 


“()- =]  (@=1), (B-1), ,F(a+2,B+2;a+B+1+ j:2) 
Á el (a+B-1), (a+B-1+ j\a+B+ j) 


Le code Mathematica suivant permet de tester numériquement ce développement : 


Clear[qPolynom2TermNO]; 
Clear[Heun2TermNOExpansion]; 
Clear[Heun2TermNOExpansionPrime]; 
Clear[Heun2TermNOExpansionPrime2]; 
Clear[Heun2TermNOExpansionED]; 


Heun2TermNOExpansion[nTerm Integer,a ,6 ,y ,z ]:= 
Sum[N[(Pochhammer[8-1,j] Pochhammer[a-1,j])/( j! Pochhammer[a+8-1,j])* 
Hypergeometric2F1[ a,6,j+a+ 6-1,2]] [j,0,nTerm]j]; 


Heun2TermNOExpansionPrime[nTerm Integer,a ,6 ,y ,z ]:= 
Sum[N[(Pochhammer[8-1,j] Pochhammer[a-1,j])/(j ! Pochhammer[a-6-1,j]) * 
(aB8Hypergeometric2F1[1+a,1+6,j+a+ 6,z])/(j-0:16-1)] [j,O,nTerm]j]; 


Heun2TermNOExpansionPrime2[nTerm Integera ,6 ,y ,z ]:- 
Sum[N[(Pochhammer[8-1,j] Pochhammer[a-1,j])/(j! Pochhammer[a 6-1, j])* 
(a(1+a)8(1+8)Hypergeometric2F1[2+ a,2+6,1+j+a+6,z]}/{(j+a+ 6- 1) (j+a+ 6))] (j,0,nTerm)]; 


Heun2TermNOExpansionPrime2[nTerm, a,6,y,z]*(y/z*2/(z -1)+(a+8-1-y)/(z -a)) Heun2TermNOExpansionPrime[nTerm, a,6,y,z]*- 
(a. *6*z-q)/(z (z -1) (z -a))*Heun2TermNOExpansion[nTerm, a, 6, y,z]; 


Prenons N-1, le polynôme de degré 2 en q dont les racines déterminent les valeurs de q est le 
suivant: 


X9- Ye, Fla Birds) à = _(r+e atj rss B+j 1) sj má 
^ €, j (y +E+j -1) Si 
le (a-1fo-248-2) — 1«& -g+2ay+(2-a\2-8) 
2ay-q*(a-248-2) — e (a-1(2-a2- p) 
1 Tores-asj)r(yse-Bsj), Tres) sei 1 (rte-a)(r+e-8), a+) 
jh IVte-a)l{+e-f) Te) s  j! (yc), e, 
Symbole de Pochhammer (y +e a),,(y+e p), (s), 


> G=l c= 


Prenons N-2, le polynôme de degré 3 en q dont les racines déterminent les valeurs de q est le 
suivant: 


q+ 

g(-10-3a(8-2)+68-a(4+67))+ 

+q(83+B(118-40)+a°(11+3B(B-4)-2a(28-5)38-4)+a(4+187 +117 )+2a(6+157-Bla+11y)+a(-4-11+ 8(2+67)))- 
(a -1f8 -1Ya -2X8 - 2a -3X8 -3) 

-2ay(? 9-8 «38 )«a( 21 «318 -? 9-38 P) 

«a^ r(I8o(y +1)+18(8 -1y +1)- ap(10+11y))- 

-6a°y(y Ay -2)-0 


jeu : ; . E 
Sc Fla piyterjz) nA USERS RIS) ERU 
j-0 Cia J (+e +J -1) $ 
1 (+e Nl -q*3 3-a|3- 
Se +e = g+a(y+2) 3«(1 ay p) S, —6€6€, -— E Nta). LUE 
2 162 


Sael ce : teca (+68) s eise 
ji (+e), 53 


Symbole de Pochhammer (y +e a), (re B), (ya e); 


à finir donner des exemples et surtout regarder attentivement la question de la convergence de 
cette série. 


Construction des équations différentielles de Heun confluentes : 


Dans ce qui suit je vais suivre les pas de l'ouvrage de S.Y.Slavyanov.W.Lay « Special Functions, 
Heun's functions » pages 97 et suivantes. 


En théorie des équations différentielles ordinaires, la confluence d'une équation s'obtient lorsque 
plusieurs points singuliers réguliers viennent à se confondre (par un processus limite par exemple) 
pour former dans l'équation résultante un (ou plusieurs) point singulier irrégulier. 


La simple confluence à partir de l'équation de Heun 


La premiére confluence que l'on applique à l'équation de Heun consiste à rapprocher le point 
singulier z-a vers z-ee. Ce dernier point z= devient alors irrégulier. On applique à l'équation de 
Heun (j'introduis la notation de Slavyanov qui comporte les rangs des points singuliers à distance 
finis, suivi du rangs du point singulier z-ee): 


vf 3.21221 1-2] (pe ET sle) 


z z-]l z-a (z-1)(z-a 
2 2 


Notation de Slavyanov Ds: o, B,7,6 ; z)=z(z 1)(z a) + (re 1)(z a)+82(e-a)+2(- a+ f+1-7-6) en pz 
z z 


Equation de Heun (a a, B,y,à ; z)»(z)) =q ylz) 


Le processus limite suivant : 


a ;4 p TL Lim -e den Lia,7,6:2) 
E E #0 € € 


zs zu y ee y: Less(: Lese if dy JIE «i4 


2 
Aa; o, y, ; z)= Lim sf LiaL7,6:2)-2(6 1) P kb (z E A cU) gas 
£20 E E dz dz 


Maarse) FEE a) PE EO 


En ce qui concerne les points singuliers réguliers z=0 et z=1, les exposants caractéristiques de 
Fröbenius restent identiques à ceux de l'équation de Heun : 


-1)+ô -1 
EME ae a ae E E 
" 5 E > : E >r=l-p >| 0 0 
o( jataz q qo=Lim z'Q(z)- Lim (z-1řQ(z)=0 1-6 


Z =) z—0 zl 1-y 


La bi-confluence à partir de la simple confluence de l'équation de Heun 


On applique le processus limite suivant à l'opérateur différentiel de simple confluence : 


2 


1a; a,y,8;z)=z(z ) 5 epe issued CD) tua 
dz dz 


1 1 : 1 1 
zO£2 a-— 82 z> Lim -ede ;10.y. £ 82) 
E & € 60 E E g 


2 
= Lim exse E D : zje: fiez-))d aes] 


£0 


= Lim zr: «( y-az 2 azl 
z 


£0 
cobi se LETE 
dz? n dz 


2 
[i3]. mm o rlir2] 1, 2,4 2 À 
L (a:a,y:z)= Lim EL (i 5: 4,7; Th sez] zs tazey) az 
d 


Wa; ay: zA) =q yle) e») 2 (Less ette 


Par une simple opération homothétique sur la coordonnée z, on parvient à la forme plus générale 
donnée par exemple par Mathematica : 


Œ = 
Due HAS) fe 2 OI d m 
RAT + zu +? RON BP TE 359 
à xy a se HI ~ a = dq 
= CE = = z 
E B g^ cg 7g 
= DL: 0282,22 6-0 


En ce qui concerne le point singulier régulier z=0, les exposants caractéristiques de Frübenius 
restent identiques à ceux de l'équation de Heun : 


pe- FEF — (n -Lim ZP()-$ 0 
me. — $a Ata >r=l-p>| 0 
o) - 82-3 q,-Lim Z'Q() 1-7 


La double-confluence à partir de la simple confluence de l'équation de Heun 


On applique le processus limite suivant à l'opérateur différentiel de simple confluence : 


2 
I?Y0 a, y,8:z)- z(z 1) ES *(y(z-1)-6 are az 
dz dz 


: 1 
z>Ź ae yy l 67 Lim pen €;0,Yy a nhi 
€ € € € 


£0 


; 4zí[z d? az az zz d Z 
-Lim 4&' —|—-l|—,—6&||Y l| +E 1 +e a 
£20 e LE dz EJE EE EE dz E 


2 2 
= Lim ke-92 z+a-ey+z(z s) mesa Eloa r +e +yzta) taz 


2 
IP?l(a : ay ; z)= Lim pex £;0,y f SE) Pate 
£0 E dz dz 


Paar Aos D En), 60 


2 
Notation a—y 78 (06,7 :2)=2 eren 2+2) +a z 
z 


ea e a ue)uo | 


dz? z dz z 


Par une simple opération homothétique sur la coordonnée z, on parvient à la forme plus générale 
donnée par exemple par Mathematica : 


a = 
| d^y(Z y ê lj)dy(Z y NS 
z-fz — Efeki) 2) E y(z)=0 

sé 2 Ub i 
y-pBy  ó-ó  £- - d-q 
B B 
d^y(Z 7 ó ~\dyZ) &zZ-q 

xi MERE) E p; 9-0 


En ce qui concerne le point singulier z=0, ce dernier n'est plus régulier. En effet à finir à recalculer 
fin du à finir 


La triple confluence à partir de la bi-confluence de l'équation de Heun 


On applique le processus limite suivant à l'opérateur différentiel de bi-confluence : 


2. 
: d 
Na; ay; z)=z ——(z-aze-y)—-az I(;z)s1 
(aia. iz)- 273 ( n ( z) 
1 2 ; 2 1 a 1 
ES TIS a p y——.t-3 => Lim e(a zQ, -5 t-z; t sz) a 2))- 


Z 
E E 
; "CNET z 1 Dies 1 
= Lim EXE re ;*€ dc E t= 
£o € €` Jdz E € E LE € 


. fz 1 \ d a(z a \d Z d? 2 d 
= Lim E 3 ste zb +EQ = + ( +a)+a z 
£0 E€ & Jdz € &'jJdz E€ dz dz 


— 


Par une combinaison d'une translation et d'une homothétie, on parvient à une forme générale de 
l'équation tri-confluente (celle donnée par Mathematica par exemple) : 


z=2À dz? dz 
y 
E Erazo 2 28 (82.2 :)»)-o 
= 8-5  poLet o go 3-12 
Y y y y Y K. 
= "(2467248 747)» (2)+(@ 2-3)»(2)-0 


à finir mettre les confluences faibles (voir Slavyanov) 


Famille des équations de Heun équivalence avec les équations dítes de Heine et Wangerin 


z " m ô E ' 1 ovz-dq m 
Equation de Heun : ? (2): fz pem zb (z)4 Lx a) t 0 


contrainte o *t*l-ycóct£— »(z) ES HeunG(a,q;0,1,7,8;z) 


1 ba 
g ENS M ad 1 UE A s 
Équation dítes de Wangerin : © Y (z)4 E eI (z): zb Ce Ad) 0 
1 1 1 1 1 ba 
QE O=T y-l ò E= o=tT= q= 
4 2 2 2 2 4 


© 
a 


= o+r+1=2=y+68+e— y(z)= HeunG Bed oT aus 
422 2 


1 ba 
z 1 L Tio 2 L 2 ' L 2 4 PN 
Équation dítes de Heine ^" ? (z^ DE Ed. EL GH z(z -1Xz - a) Gr 
OT l o l p-l1y ó -1 e=l o= t=1 q Ps 


Appendice n°1 : deux propriétés des récurrences à trois termes : 4 4, +B, d,+C,d,., =0 
- + 


Formons les deux fractions continues autour de l'indice k=k, : 


Ode 6 2 d, _ TT O d, TT eo 
Notons f, 33 fi —kk, = [I^ k « k, = [I^ 
kA k+l ko  I=ko+l ko I-k 
e m, E. 
0 k = 
dés B,+C d, B, +C fA 
k k d, 
ENE S C, e f=- C, 
^ d, dis B, + À fe 
+ B, + A, Pa k kJ k-1 
k 


Propriété n°1 : lorsque C, s'annule à l'indice ko-1 alors les coefficients du développement s'annulent 
d,=0, pour les indices k < ko. 


Démonstration de la propriété n°1 : supposons que q, , - o, alors 


k = ko -1 > A, jid 2 + Bi du + Cu ud, 


ko—l ^ ko 


Si d, ,-0-»d,,-20-5d,-0 k«k, 


=0= A 


ko da + B,, di 


Comme on suppose pour la récurrence que le coefficient qa, +o est non nul, alors la fraction 
0 

continue ft) est nécessairement nulle, même si la valeur fe, pourrait tendre vers l'infini alors : 
dd 07 07 


+ Con di 
Po ) E 0 
| B, T: 4d di, 


> dy, = 0 


La propriété n°1 est démontrée. 
Propriétés n°2 : lorsque A, s'annule lorsque k=ko+1, les coefficients d,20 s'annulent lorsque k>ko. 


Démonstration de la propriété n°2 : supposons que q, . - o, alors 


k = ko T1 A, sida, + By dia + Chotta ko2 =0= C, dua d 


uoc 0-d,-0 k»k, 


Comme on suppose pour la récurrence que le coefficient q, +0 est non nul, alors la fraction 
continue fl) est nécessairement nulle, même si la valeur ft}, pourrait tendre vers l'infini alors : 
0— ot 


Ar, +l dy +1 0 
B, a +C CE iur 
ko+l ko-HJ ko+2 ko 


I E 


La propriété n°2 est démontrée. 


L'exemple du développement en série de fonctions de Laguerre des fonctions d'ondes sphéroidales 
aplaties dans le cas oü y présente de grandes valeurs : 


4 k-—-oo 
»xG)-(-9* S a le Cra- x) Ci-?et rp) y x) 
k=—00 


Illustre bien l'une de ces propriétés. Les coefficients respectent une récurrence à trois termes de la 
forme : 


a, — (1 k uXI-k) 


B,--2 (20 7 k)- u 1) Oud, , -(A* By  x.)d, tod, —0 
X, =2U+k+u+1)(+k)+1+u 
=> À, =} C, = OG 


A, > C, 


L'annulation des termes se produit lorsque : 


Coa 2 o4, 20 2 (Ik, + u) +k) ko = 1 
Au = Aya m rk ucr2YXI-k, +2) k, = 1-2 


C 


koi 


De ces deux conditions c'est la condition =0=(1+ ko + u) + ko) ko =— qui prime de toute 


évidence et l'on peut en déduire que la sommation bilatérale devient unilatérale. 


4 k-—-oo 
»G)-(-x Ma le ya x))+ (1) 0e 70970) y x) 
k=-—I 


Appendice n°2 : deux propriétés des récurrences à quatre termes : 4 q, ,+8B, d, «C, d, +D, d, -0 
Propriété n°1 : Lorsque D, s'annule simultanément aux indices ko-1 et ko-2 alors les coefficients du 
développement d, s'annulent pour les indices k < ko si C, ,20- Si également C, ,-0 alors les 


coefficients du développement d, s'annulent pour les indices k < ko-1. 


Démonstration de la propriété n°1, premier cas : C,,«0, Supposons QUE q, , —0 er d, ,—0, 
alors : 

A, d, +B, dı +C, d, +D, d,,, =0 

k =k, -1 => A, id scs Es By, d,» sa C, dua EZ D, 4d, =0= A, adus x Ba d,» + C, dua 

Si d, , et d, ,—0-»d, ,—-0 or pourque d,—0 k<k, ilsuffit que d, ,—0 

car d, ,—d, ,—d, ,—-0-»d,.,-0 


En introduisant la fraction suivante, on peut établir les identités : 


D, 
d D at taD 
C + 
h= : Bs. : e fix : e ffa-- CT 
d, (4, fia F B.) fa t C, 4, fia + B, 4, fia + B, 
H _ . D, -l "EE "T 
Il vient pour k=ko-1 : PLI £ -0« fa #0 faa ou fes forme indéterminée 


(4. PE + B, i | = LS C, d 


et cela méme si les valeurs y, , et y, , pourraient être des formes indéterminées qui une fois 
0 0^ 


déterminées tendrait vers l'infini, zéro ou une valeur finie. En conséquence : d, ,-0- 
2 


Pour déterminer l'autre valeur en ko-2 l'identité précédente peut aussi s'écrire sous la forme : 


Df, 


fifa uc fa * Àj) 


D oft, + 


Soit pour la valeur k=ko-2 : f f, = 
ET Cath Jos + Ba 


yv te d, ,-0 


Toujours pour les mêmes raisons que f . et f, , pourraient être des formes indéterminées qui 
97 d 07 


une fois déterminée tendrait vers l'infini, zéro ou une valeur finie la valeur de 4, . s'annule. 
07 


Deuxiéme cas : C, ,-0, SUPPOSONS que d, , —o, alors : 
0 07 


A, d, ; + By d, , FC, d, +D, d,,, =0 
k =k, -1 => A, id us * Bua d, 5 + C, dua 2E D, a d,, =0= A, adus — d, 5 = 0 


k =k -2 > 4 wd WB gu BULL aD. dius de 50> dpa 50d, S0 FER el 
Et l'on démontre que 4, 


0 —2 


_0 de la même manière que précédemment. 


La propriété n°1 est donc démontrée. 


Propriétés n°2 : Lorsque A, s'annule simultanément aux indices k-kg*1 et k-kg*2, alors les 
coefficients d,=0 s'annulent lorsque k>k. 


Démonstration de la propriété n°2 : supposons que dya 
A, d; +B, dı +C, d, +D, d,,, =0 
k = k +2 => A, saa, az B,,:2 din + on 2 nz D, od ss =0= B,,:2 din al or k12 + Dy, sa uis 
SF dj 
car d,,,—d,,;—d,,,—0-»d,,,—0 


-0 et d,,,-0; alors 


et d,,,—-0-»d,,,—-0 or pourque d,—0 kk, ilsuffit que d,,,—0 


En introduisant la fraction suivante, on peut établir l'identité : 


A 


-+ +B, 
9 zd zd des * B, 4C, d, +D, Um -0e A HB +C, f e D ffa f nod 
dia da da ds ha C, * fi 4D, 
A 
g kl +B, ) 
Il vient pour k-kg*1 : f ot 


-0« f,20 f, forme indéterminée 


kotl E 


C, a + fi aDia 
et cela méme si la valeur de y, , pourrait être une forme indéterminée qui une fois déterminée 


tendrait vers l'infini, zéro ou une valeur finie. En conséquence : d, ,, -0- 
0 


Pour déterminer l'autre valeur en ko+2 l'identité précédente peut aussi s'écrire sous la forme : 


A, * B, f, 

fifa ELE 

C, + fL D, 

Aa a B, fia I An E d, 


Soit pour la valeur k-ko*2 : f, f, =- 
v 3 C a t ÁaDa C, ua * fai d, 


Toujours pour les mêmes raisons que f . pourrait être une forme indéterminée qui une fois 


déterminée tendraient vers l'infini, zéro ou une valeur finie dés lors la valeur de a, . s'annule. 
0 


La propriété n°2 est donc démontrée. 


Appendice n*3 : deux propriétés des récurrences à cinq termes : A, d, , * B, d, 4 C, d, 4 D, d, +E, d, -0 
k^ k- - + + 


Propriété n°1 : Lorsque D, et E;s'annulent simultanément à l'indice ko-1 et E,s'annule à l'indice ko-2 
alors les coefficients du développement s'annulent d,-0, pour les indices k < ko. 


Démonstration de la propriété n°1 : supposons que d, | —-0 er d, ,-0, alors 


A,d, +B, dy +C, dy +D, dya +E, d,,, =0 

k-ky—1— Ay -1dr -3 + Bp-i d, 5 + C, id a t D, dy, * E, ada =0= Ay ad, + Bp-i d, -2 + Chide 
Si d, ,-0 et dp- =0— d, ,-3=0 

Pour que d,=0 k<k, ilsuffit que d, ,—0 

k =k, -2 > Ar -2d k -4 * B, -2 d, 5 + C, ad, -2 t D, ad, a t E, -2dr = 0 = A, cada t E, ad, 

Comme E, ,-0-»d, ,-0-»d,-0 k<k, 


En introduisant la fraction suivante, on peut établir l'identité : 


d Pi E 
d = k a = k+l 
di, À, fiafia y B, Ja + C, 
D, ,+ Er 
Il vient pour k=ko-1 : ME: 


ni =0 f, #0 fan ou f,, forme indéterminée 


À: fihi * B, , faa t Cpa 


et cela méme si les valeurs f , et y, , pourraient être des formes indéterminées qui une fois 
+ x 0— 


déterminées tendrait vers l'infini, zéro ou une valeur finie. En conséquence : d, ,-0- 
z 


Pour déterminer l'autre valeur en ko-2 l'identité précédente peut aussi s'écrire sous la forme : 


D, fin +E, 
C, + fia (4, fa + B,) 


Sefin = 


Soit pour la valeur k=ko-2 : pem Dfii Fus E dus 


C, a + fas 4... fa T B, | i Cp- * halia ne + 85] d, 


Toujours pour les mêmes raisons que y, . et f, , pourraient être des formes indéterminées qui 
0— * ; ES 


une fois déterminée tendrait vers l'infini, zéro ou une valeur finie la valeur de 4, . s'annule. 
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La propriété n°1 est donc démontrée. 


Propriétés n*2 : Lorsque A, et B, s'annulent simultanément lorsque k-ko*1, et que A, s'annule 
lorsque k=k,+2, alors les coefficients d,=0 s'annulent lorsque k>ko. 


Démonstration de la propriété n°2 : supposons que d. 


150: et pm 0x alors 
En d, 5 +B, d, , +C, d, +D, di +E, di, = 0 


k = ko +1 => A, ada + Ba d, + C, ad ua + D, di, + E, adus =0= C, ad na ze D, adio + Ex adu 
Si d, a =0 et d, i5 =0—d,,,: = 0 
Pour que d;=0 k>k, ilsuffit que d,+4=0 
k = ko +2 => A, adl, * B,,:2 dia + o kot2 * D, aad qs H Ead usa =0= A, sad, * Es sadysa 
Comme 4,,,=0—=4d,,,=0—=d,;=0 kk, 
En introduisant la fraction suivante, on peut établir l'identité : 
À 

NP. ot 2 4 ue 
h= : 4, i B, 4 C, : t D, Hy k H 20e —-4 B, +C, f + D, f fia E, fff > f us 

d, dya da d, , d, , fa C, + fin (D, +E, faa) 


Il vient pour k=ko+1 : 


p -0cf,*0  f,, ou f, formeindéterminée 
| Cu + ne (D... + Epa fus) : ' i 


et cela méme si les valeurs f , et f,., pourraient être des formes indéterminées qui une fois 
0 0 


déterminée tendrait vers l'infini, zéro ou une valeur finie. En conséquence : d, ,, -0- 
0 


Pour déterminer l'autre valeur en ko*2 l'identité précédente peut aussi s'écrire sous la forme : 


À, +B, fi: 
C, + f (D, +E, fan) 


fifin T 


P A, ., t B, A 
Soit pour la valeur k=ko+2 : faha h2 T "Y N ho E d, ia 


Cu fida *E, 2 fua) C, Tas DE *E, 2 dal d, 


Toujours pour les mêmes raisons que f . et pf, , pourraient être des formes indéterminées qui 
A 0 M 0 
une fois déterminées tendraient vers l'infini, zéro ou une valeur finie dés lors la valeur de q, , 
0 


s'annule. 


La propriété n°2 est donc démontrée. 


